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I.  INTRODUCTION 

Seals  develop  fluid-struoturo-lnteraotion  forces  due  to  the  relative 
notion  of  the  rotor  and  housing  at  seal  locations,  and  these  forces  have 
been  shown  to  have  a  significant  influence  on  the  rotordynamio  stability 
and  response  characteristics  of  high  performance  turbomachinery.  The 
present  researoh  effort  was  initiated  in  September  of  1982,  and  is  a 
oombined  computational-experimental  program  to  develop  experimentally- 
validated  computational  approaohes  for  predicting  seal  forces.  More 
speoif ioally,  computational  approaches  are  being  developed  to  define  the 
force  coefficients  for  the  foroe-motion  relationship: 
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where  (X,Y)  are  components  of  the  relative  displacement  veotor  between  the 
rotor  and  the  stator  at  the  seal  location,  and  (Fy,  Fy)  ***•  components  of 
the  reaction  force  resulting  from  transient  pressure  and  velooity 
distributions  within  the  seal.  The  computational  approaohes  which  are 
being  developed  in  this  study  will  define  the  pressure  and  velocity 
distributions  within  the  seal  whioh  result  from  relative  rotor-stator 
motion.  The  computational  approaohes  are  to  be  validated  by  comparison  to 
experimental  measurements  of  the  following  quantities: 

(a)  transient  measurements  of  the  displacement  (X,  Y)  and  foroe 
(Fx,  py )  components,  whioh  results  in  calculated  values  for  the 
rotordynamio  (stiffness  and  damping)  coefficients  of  Eq.  (1),  and 

(b)  laser  velooimeter  measurements  of  the  velocity  distributions 
within  labyrinth  seal  oavities. 

Dr.  Childs  and  Dr.  Rhode  are  responsible,  respectively,  for  the 


experimental  and  computational  developments.  In  the  following  two 
chapters,  reviews  are  provided  of  the  research  aooomplishments  with  respeot 
to  these  two  complementary  faoets  of  the  research  projects. 

II.  EXPERIMENTAL  DEVELOPMENTS 

A,  Introduction 

The  experimental  developments  of  this  program  are  Jointly  supported  by 
NASA  and  AFOSR.  NASA  funds  ($323,749;  1  July  1981  to  31  Deoember  1984) 
have  been  used  to  develop  the  test  faoility  and  test  apparatus  and  to 
support  analytloal-computatlonal  developments.  NASA  officials  are 
primarily  interested  in  annular  seals  having  oonstant-olearanoe  or 
oonvergent-taper  geometries  with  smooth  and  honeyoomb  stator  surfaoes. 
Their  interest  in  this  type  of  seal  arises  from  applications  in  the  turbine 
areas  of  the  high-pressure  turbopumps  of  the  Spaoe  Shuttle  Main  Engine. 
This  status  report  reviews  the  progress  of  the  complete  program  without 
distinguishing  between  souroes  of  support. 

B.  Research  Program  2  Review 

During  the  past  year,  the  test  program  has  )een  in  a  "shakedown"  mode 
involving  the  following  developments: 

(a)  Aooelerometers  were  installed  on  the  stator  to  oorreot  for 
aooeleration-induoed  forces  whioh  are  measured  by  the  transient  load 
cells.  Software  and  eleotronios  have  been  developed  to  implement  this 
correction. 

(b)  One  of  the  hollow  roller  bearings  whioh  support  the  test  shaft  in 
its  housing  was  installed  improperly  by  the  original  machine  shop  and 
failed.  Replacement  units  were  installed  properly  and  have  worked 
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satisfactorily. 

(o)  The  original  coolant  and  lubricant  loop  for  the  radial  bearings 


and  the  thrust  bearing  was  found  to  be  inadequate  and  was  replaced. 

(d)  A  stepper  motor  was  added  to  provide  remote  external  control  for 
vertical  positioning  of  the  test  section  seal. 

(e)  A  4:1  belt-pulley  speed  inoreaser  was  added  to  the  eleotrio-motor 
drive  system  to  provide  a  theoretical  top  speed  of  21,500  opm.  This 
eliminated  the  need  for  the  installation  of  an  air  turbine  whioh  was 
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originally  provided  for  higher  speeds.  The  belt-drive  system  also 
eliminated  the  undesirable  restraint  of  the  coupling  which  was 
originally  provided  as  a  low-speed  drive  system.  The  pulley  system  has 
a  marginal  but  adequate  lifetime  out  to  20,000  rpm  when  transmitting 
up  to  30  hp. 

(f)  Extensive  software  has  been  developed  to  expedite  the  execution  of 
testing  and  the  processing  and  presentation  of  test  data.  As  a  result, 
150  data  points  oan  generally  be  secured  in  under  three  hours. 

(g)  An  approximate  analysis  was  developed  [1]  based  on  Iwatsubo’s  work 
[2]  for  the  prediction  of  leakage  and  the  rotordynamio  coefficients  of 
Eq.  (1)  for  labyrinth  seals.  A  oopy  of  reference  C 1  3  is  attached.  In 
comparison  to  previously  published  data  of  Waohter  and  Benokert  [3], 
the  analysis  provides  satisfactory  results.  However,  Waohter  and 
Benchert  provide  only  direot  and  cross-coupled  stiffness  data  (K  and  k 
of  Eq.  (1)),  with  no  damping  measurements.  Further,  the  data  only 
applies  for  fluid  prerotation  without  shaft  rotation,  or  shaft  rotation 
without  fluid  rotation.  As  specified  in  [1],  the  approximate  analysis 
is  in  rough  agreement  with  teeth-on-stator  data  but  is  completely  at 
odds  with  the  published  teeth-on-rotor  results. 
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(h)  The  following  seal  configurations  were  tested: 

(i)  smooth,  oonstant-olearanoe, 

(li)  smooth,  oonvergent-tapered, 

(iii)  smooth-rotor,  honeycomb  stator, 

(iv)  constant-clearanoe  labyrinth  with  teeth-on-stator,  and 

(v)  oonstant-olearanoe  labyrinth  with  teeth-on-rotor . 

The  honeyoomb  seal  was  provided  by  NASA  and  Rocketdyne  and  is  the 
turbine  interstage  seal  for  the  HPOTP  (High  Pressure  Oxygen  Turbopump) 
of  the  SSME  (Spaoe  Shuttle  Main  Engine).  The  two  labyrinth  seals  are 
illustrated  in  figure  1  and  use  cavity  and  tooth  dimensions  based  on 
0.  E.  aircraft-gas-turbine  design  practice. 


C.  Labyrinth  Seal  Teat  Results 

The  test  apparatus  and  faoility  are  described  in  detail  in  Nicks' 
thesis  [4],  which  is  attaohed  as  part  of  this  report.  The  test  apparatus 
has  the  capability  to  identify  the  seal  rotordynamio  coefficients  of  the 
general  linearized  model 
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for  small  motion  about  an  arbitrary  static  eccentricity  position  defined  by 

®o,  e.g.,  Xo  ■  eo,  Y o  ■  This  is  in  oontrast  to  the  model  of  Eq.  (1) 
whioh  applies  for  small  motion  about  a  oentered  position. 

The  following  parameters  oan  be  controlled  in  testing: 

(a)  Pressure  ratio.  The  seal  discharge  pressure  is  atmospheric,  and 
the  supply  pressure  oan  be  increased  from  zero  to  approximately  90  psi. 

(b)  Running  speed.  The  running  speed  range  for  the  present  test 
apparatus  is  zero  to  8000  rpm. 
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(o)  In3 et  Tangential  Velocity.  Swirl  rings  are  provided  at  the  seal 
Inlet  to  prerotate  the  flow  either  in  the  same  direction  or  an  opposite 
direction  to  rotation.  Also,  fluid  oan  be  introduced  with  no 
intentional  prerotation. 

All  testing  is  done  with  air  and  includes  static  measurements  of  the 
entrance  and  exit  pressures  and  temperatures,  and  pressure  measurements  in 
each  labyrinth  cavity.  The  testing  reported  here  _is  for  small  motion  about 
a  oentered  position. 

Figure  2  illustrates  the  test  results  for  a  tooth-on-rotor  labyrinth 
seal  configuration  with  fluid  prerotation  in  the  direction  of  shaft 
rotation.  The  independent  variable  is  running  speed  and  families  of  curves 
are  presented  for  various  pressure  ratios.  Observe  that  all  of  the 
coefficients,  exoept  for  the  oross-ooupled  damping  coefficients,  are 
relatively  insensitive  to  running  speed.  The  oross-ooupled  damping 
coefficient  displays  the  most  erratio  behavior  with  respeot  to  its 
dependency  on  running  speed  and  pressure  ratio.  Again,  from  figure  2, 
observe  the  general  monotonio  inoreaae  of  the  rotor  dynamic  ooeffioients 
with  Increasing  pressure  ratios. 

From  a  rotordynamios  viewpoint,  the  direct  damping  and  oross-ooupled 
stiffness  ooeffioients  are  quite  significant  in  determining  the  stability 
and  peak  response  characteristics  of  rotors,  The  direot  stiffness  and 
oross-ooupled  damping  terms  have  a  minute  influence  on  oritioal  speed 
locations  and  are  of  markedly  less  interest.  Figures  three  and  four 
illustrate  the  dependence  of  these  coefficients  on  pressure  ratio.  The 


three  curves  on  eaoh  of  these  figures  are  for  no  fluid  prerotation  and 
prerotation  in  the  same  and  opposite  direction,  respectively,  to  shaft 
rotation.  The  curves  oonneot  the  average  value,  with  respeot  to  speed,  for 


Direct  stiffness  results  for  teeth- on  rota*  labyrinth  seal 

configuration  with  fluid  prerotation  in  the  direction  of  shaft 
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FIGURE  3(b).  Direct  daaping  versus  pressure  ratio  far  three  inlet-tangential 
velocity  conditions;  teetb-on-rotor  configuration. 


FIGURE  4(b).  Direct  damping  versus  pressure  ratio  for  three  inlet  tangential 
velocity  conditions;  teetb-oh-stator  configuration. 
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each  coefficient.  The  horizontal  bars,  whioh  oonnect  the  vertical  lines 
through  eaoh  data  point  show  the  extreme  values  due  to  speed  dependency. 

With  respect  to  the  oross-coupled  stiffness,  the  results  of  figures  3 
and  4  show  a  strong  influence  on  fluid  prerotation  with  k  increasing  and 
decreasing  due  to  fluid  prerotation,  respectively,  in  the  direotion  of 
shaft  rotation  and  opposed  to  shaft  rotation.  Comparatively  small  values 
result  for  no  shaft  rotation.  Cross-coupled  stiffness  values  are  generally 
smaller  for  teeth-on-stator  than  teeth-on-rotor  configurations. 

The  direct  damping  value  is  loss  sensitive  to  fluid  prerotation.  It 
is  generally  largest  for  fluid  prerotation  against  shaft 
rotation,  moderately  smaller  for  prerotation  in  the  some  direotion  as 
rotation,  and  smallest  for  no  fluid  prerotation.  Observe  that  direot 
damping  increases  with  the  pressure  ratio  and  would  have  an  approximate, 
extrapolated  zero  value  for  a  pressure  ratio  of  one. 

The  test  results  provided  in  figures  2  through  4  are  vastly  better 
than  any  whioh  have  previously  been  published.  They  have  only  recently 
been  obtained  and  will  obviously  require  substantial  additional  study.  The 
correlation  with  the  analysis  of  reference  [1],  developed  in  this  study,  is 
disappointing.  Speoifloally,  measured  oross-ooupled  stiffness  coefficients 
are  muoh  less  sensitive  to  running  speod  than  predicted,  and  measured 
damping  values  are  on  the  order  of  four  times  higher  than  predictions.  As 
disouased  in  the  following  seotlon,  a  "new  and  improved"  but  still  quite 
simplified  solution  will  be  developed  in  the  ooming  year  whioh  should 
substantially  improve  the  correlation  between  theory  and  experiment. 


D.  Hssearoh  Program  2  PI*1**  for  Current  Tear 

The  following  specif io  tasks  will  be  completed  in  current  contraot 
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time  period: 

(a)  Teat  Apparatus  Development.  The  labyrinth  seals  leak  a  great  deal 
leas  than  the  annular  seals  which  have  been  tested  in  the  past  for 
NASA.  As  a  oonaequenoe  the  vortex  flow  meter  whioh  was  originally 
purohased  has  proven  to  be  oversizod  and  will  shortly  be  replaoed  by  an 
appropriately  sized  turbine  flow  meter.  In  addition,  the  inlet  swirl 
rings  were  also  oversized  to  achieve  sufficiently  high  tangential 
velocities  and  are  hei'nj.  supplemented  with  new  units.  Finally,  various 
redesigns  are  underway  to  yield  a  sharp  lnoreaae  in  running  speed, 
probably  to  something  on  the  order  of  15.000-16,000  opm. 

(b)  Computational  Developments.  As  noted  earlier,  a  simplified 
analysis  has  been  developed  for  the  prediction  of  rotordynamio 
coefficients  of  labyrinth  seals  [1],  and  has  not  proven  to  be 
particularly  satisfactory.  The  motivation  for  developing  this  type  of 
analysis,  as  oompared  to  the  oomplete  three-dimensional  Navler-Stokes 
solution  of  Dr.  Rhodes,  is  quite  simply  oost  and  oonvlence.  A 
relatively  simplified  analysis  i»  needed  for  design  purposes.  A  ttnew 
and  improved"  simplified  analysis  is  currwntly  under  development  based 
on  the  work  of  Jenny,  Wyssman,  and  Pham  [5].  We  propose  to  "calibrate" 
this  model  by  comparison  to  Dr.  Rhodes'  more  comprehensive  solution. 

(o)  hSA  Developments.  Dr.  Morrison  and  his  graduate  student,  Mr.  Mike 
Long,  will  oontlnue  to  work  at  achieving  velooity  measurements  within  a 
labyrinth  oavity. 

B.  GRADUATE  STUDENT  DEVELOPMENT 

Mr.  Joseph  Soharrer  has  been  working  on  this  project  since  1  January 

1983.  Mr.  Soharrer  concluded  his  B,  S.  degree  in  Mechanloal  Engineering  at 
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Northern  Arizona  University  in  July  of  1982  (in  three  years).  He  has  been 
tooepted  into  the  AFRAPT  program  by  TAMU  and  spent  the  summer  of  19811  with 
General  Eleotrio  at  Lynn,  Massachusetts.  Mr.  Soharrer  has  been  responsible 
for  a  great  deal  of  the  experimental  development  work  on  this  projeot,  and 
has  reoently  assisted  in  the  development  of  an  Improved  prediction  model 
for  rotordynamio-ooefficienta  of  a  labyrinth  seal  [1].  Mr.  Mike  Long  has 
been  working  with  this  projeot  for  the  past  year  with  Dr.  Gerald  Morrison 
on  the  development  of  LDA  measurement  techniques  for  labyrinth  seals.  He 
spent  the  summer  working  with  Garrett  at  Phoenix,  Arizona. 


HI.  COMPUTATIONAL  DEVELOPMENTS 

A.  Introduction 

Good  progress  has  bean  mads  In  developing  an  advanced,  user-simplis¬ 
tic,  experimentally-validated,  finite  difference  computer  code  for 
predicting  the  pressure  and  velocity  fields  and  subsequently  the  four 
rotordynamic  force  coefficients  of  Eq.  (1).  The  problem  of  subsonic 
air  flow  through  a  labyrinth  seal  is  clearly  elliptic,  as  recirculation 
tones  have  been  observed  in  flow  visualization  experiments  [6].  Also, 
it  is  three-dimensional  because  of  rotor  eccentricity  effects  which 
are  of  primary  concern.  It  has  been  assumed  that  the  rotor  Instability 
is  manifested  as  a  circular  precossional  orbit  of  frequency  w  about 
the  seal  canter  position.  The  desired  solution  will  be  obtained  from 
the  steady-state  form  of  the  governing  flow  equations  because  a  steady- 
state  flowflald  processes  along  with  the  rotor  in  a  r-0  plane. 

Once  the  predicted  pressure  distribution  (for  a  given  orbital 
velocity  u)  acting  on  the  rotor  is  obtained,  it  will  be  employed  to 
calculate  the  resulting  and  F0  forces.  A  second  computation  of  the 
pressure  and  rotor  force  components  at  a  different  value  of  u>  finally 
allows  each  of  the  all-important  stiffness  (K,k)  and  damping  (C,c)  coef¬ 
ficients  to  bo  calculated  as  needed  for  dynamics  calculations  of  the 
rotor. 

B.  False  Diffusion 

Until  recently,  the  upwind  differencing  scheme  has  been  the  most 
effective  method  for  obtaining  a  stable  numerical  solution  of 
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convection-dominated  flows.  However,  this  scheme  la  only  first-order 
accurate  and  may  Introduce  serious  truncation  errors  known  as  false 
diffusion.  This  is  especially  true  for  flows  where  convection  dominates 
(l.e. •  where  the  grid  Peclet  number  |Pe|=|v|A/r^  exceeds  2.0)  with  sub¬ 
stantial  streamline-to-grid  skewness  and  diffusive  transport  normal  to 
the  flow  direction. 

False  diffusion  results  in  an  overly  diffusive  solution.  A  turbu¬ 
lent  recirculating  flow,  such  as  that  within  some  labyrinth  seals,  is 
particularly  susceptible  to  the  effects  of  false  diffusion  because  of 
its  regions  of  considerable  velocity  gradient  and  streamlino-to-grid 
skewness.  One  way  to  prevent  false  diffusion  is  to  utilize  a  sufficiently 
fine  grid  so  that  central  differencing  for  the  convective  terms  may  be 
used  throughout  the  calculation  domain.  However,  this  is  quite  expensive 
for  two-dimensional  flows  and  prohibitive  for  three-dimensional  ones. 

The  Hybrid  upwind/central  differencing  scheme  was  originally 
embodied  in  the  incompressible  flow  code.  Using  this  method,  the  con¬ 
vective  terms  are  formulated  by  upwind  differencing  for  a  grid  Peclet 
number  larger  than  2.0;  otherwise,  asntral  differencing  is  used.  This 
differencing  scheme  has  yielded  solutions  in  which  false  diffusion  has 
entirely  obscured  physical  diffusion  [7,8]. 

The  recent  QUICK  scheme  derived  by  Leonard  [9]  eliminates  false 
diffusion.  This  scheme  has  been  implemented  in  an  original  and  highly 
stabilizing  way.  It  eliminates  false  diffusion  by  incorporating  a  three- 
point  upwind-shifted  quadratic  interpolation  formula.  Rhode,  et  al,  [10] 
showed  comparisons  between  the  Hybrid  and  QUICK  differencing  schemes  for 
lncomprtssible  flow  in  labyrinth  ssul  oavltlss.  It  waa  shown  that  the 
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QUICK  scheme  yields  a  grid-independent  solution  using  a  considerably 
more  coarse  grid  than  does  the  Hybrid  scheme,  thereby  considerably 
reducing  the  computational  effort  required.  In  that  concentric-rotor 
eeal  flow  problem,  it  was  found  that  a  solution  of  given  accuracy  was 
obtained  using  QUICK  while  consuming  only  44%  of  the  execution  time 
required  by  the  state-of-the-art  Hybrid  scheme.  This  advantage  will 
certainly  be  more  dramatic  for  the  three-dimensional  (eccentric-rotor) 
seal  flow  problems  which  lie  ahead.  This  computational  advancement  was 
also  presented  by  Rhode,  et  al.  [11]  for  corresponding  compressible 
flows. 

Both  the  incompressible  and  the  compressible  flow  labyrinth  seal 
papers  above  Included  the  complete  distribution  of  numerous  quantities 
revealing  the  detailed  flow  character  throughout  a  realistic  labyrinth 
seal  cavity.  Previously  these  detailed  results  were  largely  unavailable. 
Those  results  serve  to  enhance  one's  fundamental  understanding  of  the 
compressible  flow  field  in  a  labyrinth  seal  which  generates  the 
rotordynamic  instability  forces  of  interest  to  gas  turbine  manufacturers. 

C.  Initial  Compressible  Flow  Prediction  Test  Case  of  Concentric-Rotor 
(2-D)  Computer  Code 

Due  to  the  late  arrival  of  the  Laser-Doppler  anemometer,  detailed 
measurements  of  compressible  flow  in  a  labyrinth  seal  are  yet  to  be  ob¬ 
tained.  Thus,  a  brief  prediction  assessment  using  a  compressible 
recirculating  flow  problem  was  conducted  before  developing  the  final 
three-dimensional  version  of  the  code.  A  short  literature  search  led 
to  the  selection  of  the  sudden- expansion  dump  combustor  flowfield 


measured  by  Drewry  [12].  He  measured  stagnation  and  static  pressure 
as  well  as  static  temperature  at  five  axial  stations. 

The  test  facility  employed  by  Drewry  is  shown  schematically  in 
Fig.  5.  High  pressure  air  enters  a  large  plenum  chamber  which  exits 
into  the  63.5  mm-diameter  inlet  duct.  After  passing  through  this  317.5 
mm-long  duct,  it  enters  the  combustor  test  section  which  exhibits  a 
diameter  and  length  of  97.5  mm  and  381  mm,  respectively.  Finally,  the 
flow  discharges  through  a  nozzle  of  63.5  mm-diameter  into  the  exhaust 
duct. 

This  flowfleld  is  characterized  by  both  the  Reynolds  number  and  the 
Mach  number.  Based  on  inlet  duct  conditions,  these  parameters  exhibit 
values  of  1.42x10^  and  0.67,  respectively.  The  inlet  static  pressure 
and  temperature  are  156  kFa  (absolute)  and  249. 8°K,  respectively, 

The  inlet  flow  measurements  are  perhaps  most  accurate  at  the 
measurement  station  approximately  105.0  mm  upstream  of  the  duct  expan¬ 
sion  (combustor  inlet);  hence,  this  location  was  chosen  for  the  upstream 
computational  boundary  where  a  uniform  profile  is  assumed  for  the  axial 
velocity  and  temperature.  Values  of  x  in  the  figures  are  measured  from 
this  reference  location.  The  domain  extends  downstream  348.0  mm  from 
the  expansion.  Solutions  were  obtained  using  the  QUICK  scheme  with  a 
34  x  30  grid  as  well  as  the  Hybrid  scheme  with  a  53  x  30  grid. 

The  overall  flow  pattern  is  indicated  by  the  predicted  streamline 
plot  shown  in  Fig.  6.  This  result  was  obtained  using  the  QUICK  scheme 
with  the  34  x  30  grid.  The  stream  function  ie  nondimensionalized  by 
the  inlet  mass  flow  rate.  The  recirculating  mass  flow  rate  within  the 
recirculation  zone  is  approximately  10%  of  the  inlet  value.  Drewry' s  oil 
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surface  flow  visualization  results  indicated  that  the  dividing  stream¬ 
line  (tfi*-  1.0  )  reattachment  occurs  at  x/D-1,99.  This  length  is 
under-predicted  here  by  only  2.2%.  The  53  x30  Hybrid  scheme  solution 
(not  shown)  is  almost  identical,  except  in  that  case  the  reattachment 
length  is  tinder-predicted  by  6.2%. 

The  radial  profile  plots  presented  in  Figures  7  through  10  reveal 
the  degree  of  agreement  with  measurements.  The  QUICK  differencing 
scheme  was  used  along  with  the  34  x  30  grid  in  obtaining  these  predic¬ 
tions.  Figure  7  shows  predicted  and  measured  time-mean  axial  velocity 
profiles,  nondimensionalized  by  the  spatially-averaged  inlet  duct  velocity 
U*212.3  m/s.  Observe  the  excellent  agreement  at  the  first  and  third 
axial  stations  in  the  presence  of  significantly  lass  agreement  at 
the  second  station.  This  is  somewhat  surprising  as  the  large  3u/8r 
gradient  at  the  first  station  appears  to  represent  the  most  substantial 
prediction  challenge.  The  discrepancy  at  the  second  station  is  attri¬ 
buted  to  measuremiMt  inaccuracy  due  to  the  use  of  non-dlractional 
pressure  probes.  This  occurs  for  r/D>0.3  where  the  streamlines  were 
seen  to  exhibit  non-zero  slope;  hence,  relative  to  the  local  flow 
direction,  the  probe  was  oriented  at  a  significant  pitch  angle  rather 
than  facing  directly  into  the  flow  as  desired. 

Static  pressure  is  given  in  Fig,  8.  Each  radial  profile  is  nearly 
uulform,  and  the  agreement  with  measurements  is  excellent.  The  dis¬ 
crepancy  is  lass  than  3.0%  almost  everywhere.  Figure  9  shows  profiles 
of  static  temperature.  The  slight  nonuniformity  seen  in  these  profiles 
reflects  the  large  velocity  change  occurring  over  the  free  shear  layer 
thickness  in  this  flow  of  essentially  constant  stagnation  enthalpy. 
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Excellent  agreement  with  measurements  1b  again  obtalnad. 

Pradlctad  density  profila*  also  exhibit  this  nonuniformity  as 
shown  In  Fig.  10.  Due  to  tha  lnvaraa  ralationahlp  between  density  and  f\ 
temperature,  tha  radial  gradient  la  opposite  to  that  of  temperature, 

D.  Prediction  of  Compressible  Flow  In  a  Concentric-Rotor  Labyrinth  Seal 

The  convergence  rate  and  other  operational  characteristics  of  tha 
two-dimensional  (concentric-rotor)  compressible  flow  code  ware  evaluated 
via  prediction  of  a  generic  labyrinth  seal  of  the  same  geometry  which 
was  computed  earlier  for  Incompressible  flow,  This  also  serves  to 
enhance  one's  Insight  Into  tha  compressible  flow,  eccentric-rotor  flow 
field  which  produces  the  all-important  rotordynamic  forces.  Tha  configu¬ 
ration  of  the  labyrinth  seal  considered  Is  depicted  in  Fig.  11. 

The  seal  configuration  exhibits  a  cavity  axial  length  L-1.113  mm, 
stator  wall  radius  Ri"42.89  mm,  tooth  radial  clearance  c«0.216  mm,  and 
radial  distance  from  cavity  base  to  stator  wall  d-1,105  mm.  The  shaft 

speed  ft"35,410  rpm,  the  mean  cavity  inlet  temperature  T"294.4*K,  and 

2 

the  cavity  Inlet  static  pressure  on  the  stator  wall  Pow»3. 06x10  IcPa 
(absolute).  Results  are  presented  for  the  case  in  which  the  cavity 
inlet  Mach  number  M^-0.65  using  the  tooth-clearance  bulk  axial  velocity. 
Other  dimensionless  parameters  are  Re 32Uc/v-2. 07x10^  and 

L  3 

TaS(Wd/v)(d/rih)  -4.16x10  .  Solutions  were  also  obtained  for  inlet  Mach 
number  values  of  0.2,  0.5,  and  0,72  using  the  same  inlet  pressure,  den¬ 
sity,  and  temperature.  Comparisons  between  the  four  solutions  are 
discussed  subsequently. 

Measured  distributions  of  compressible  flow  variables  within  a 
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labyrinth  iuI  cavity  ara  not  yat  available.  For  tha  pravloua  inoom- 
praaalblo  flow  raaulta  [10],  thaaa  Inlet  profiles  vara  calculated  along 
with  tha  solution  by  setting  tha  inlet  value  of  each  variable  (except 
pressure)  equal  to  tha  corresponding  outlet  value  at  each  numerical 
iteration.  This  was  based  on  the  assumption  that  the  incompressible 
flow  is  streetwise  periodic,  However,  this  assumption  Is  not  valid  for 
tha  present  variable-density  flow  situations,  For  the  present  case, 
tha  dimensionless  inlet  profiles  of  the  incompressible  flow  solution 
were  adopted  as  the  dimensionless  inlet  boundary  values  for  compressible 
flow  through  tha  identical  seal  cavity, 

Figures  12  through  17  exhibit  predictions  using  a  33  x  31  grid 
employing  the  QUICK  scheme,  This  arrangement  was  previously  shown  [10] 
to  yield  accurate  solutions  of  the  difference  equations  for  the  incompres¬ 
sible  flow  case.  The  convergence  criteria  for  this  study  ranged  as  low 
as  0.1X  for  the  normalised  sum  of  residual  source  magnitudes. 

Figure  12  shows  the  predicted  dimensionless  streamline  pattern. 

The  stream  function  is  nondimensionallsed  by  the  leakage  mass  flow  rata, 
The  free  sheer  layer  emanating  from  tha  separation  comer  gives  rise  to 
the  expected  recirculation  sona  in  the  cavity,  There  is  no  indication 
from  predictions  that  an  additional  small  recirculation  sons  exists  in 
the  lower  comers  of  the  cavity  near  the  base  of  each  tooth. 

Although  atraamllne  patterns  for  other  oases  ara  nearly  identical, 
an  interesting  observation  can  be  made  from  comparing  Fig.  12  with  stream¬ 
line  plots  for  M^-Q.2,  0,3,  and  0,72,  The  recirculation  strength  as  a 
percentage  of  leakage  mass  flow  rata  decreases  slightly  for  incrsasing 
,  This  rslative  strength  for  the  M^O.2  case  Is  40. 2%  and  it  ataadily 
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decreases  to  35.4%  for  M^-0. 72.  For  comparison,  the  flow  of  liquid 
hydrogen  through  the  Identical  seal  cavity  at  Ra^-lxlO^  and  Ta*l .3x10^ 
produces  a  recirculation  strength  of  39.0%. 

The  predicted  radial  profile  plots  pertain  to  the  M^-0.65  case  and 
indicate  numerous  previously  unavailable  details.  Corresponding  results 
for  other  Mach  number  cases  are  vary  similar.  One  exception  is  a  some¬ 
what  more  pronounced  variation  in  axial  velocity  and  pressure  for  the 
higher  cases  in  the  region  above  the  downstream  tooth.  Figure  13 
shows  axial  velocity  nondimenaionalised  by  the  cavity  inlet  bulk  value 
U«223.5  m/s.  The  solution  indicates  the  expected  large  value  of  3u/Sr 
near  the  separation  corner.  This  reveals  a  high  level  of  turbulence 
energy  generation  in  this  thin  shear  layer  region. 

For  the  higher  cases,  the  large  pressure  decrease  from  inlet 
to  outlet  of  the  cavity  yields  a  corresponding  density  decrease.  This, 
in  turn,  results  in  a  convective  acceleration  of  the  leakage  flow  in 
accordance  with  mass  conservation.  This  effect  is  shown  in  Fig.  14  for 
each  M^.  As  expected,  this  Mach  number  effect  sharply  Increases  with 
increasing  M^, 

The  swirl  velocity  component  is  shown  in  Fig.  15  in  which  the  shaft 
circumferential  velocity  is  the  nondimensionalizing  factor.  The  distri¬ 
bution  is  nearly  uniform  at  a  value  of  0.72.  The  corresponding  liquid 
hydrogen  case  [10]  exhibited  uniform  profiles  with  a  value  of  0.65. 

Not  shown  here  are  very  sharp  gradients  near  each  surface. 

Figura  16  shows  the  preasura  distribution  relativs  to  the  inlet 
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stator  well  preasura  F  ,  which  for  this  case  is  3.06x10  Pa  (absolute). 

ow 

It  la  normalized  by  the  cavity  inlet  dynamic  pressure.  The  pressure 


decreases  slowly  until  x/L-0. 5  and  then  increases  slightly  until  Approxi¬ 
mately  x/L"0.75,  Finally,  it  decreases  sharply  as  the  flow  accelerates 
over  the  downstream  tooth.  The  sharp  pressure  peak  observed  at  x/l>0.85 
results  from  flow  stagnation  on  the  downstream  tooth.  For  the  given 

leakage  flow  rate,  the  overall  bulk  pressure  drop  across  the  cavity  is 

a 

approximately  3.5x10  Fa  (0,34  atm). 

Dimensionless  turbulence  kinetic  energy  is  plotted  in  Fig.  17.  As 
mentioned  previously,  it  is  the  large  value  of  &u/dr  In  the  free  shear 
layer  which  produces  intense  turbulence  energy  in  that  region,  This 
promotes  the  required  large  bulk  pressure  drop.  As  expected,  the  greatest 
turbulence  energy  value  occurs  near  the  stagnation  point  on  the  downstream 
tooth,  The  recirculation  region  effectively  acts  as  a  turbulence  energy 
sink  as  the  figure  indicates  a  lack  of  turbulence  generation  there. 

E.  Eccentric-Rotor  (3-D)  Computer  Code  Development 

The  final  major  task  in  developing  the  required  eccentric-rotor 
computer  code  entails  the  extension  of  the  concentric-rotor  compressible- 
flow  code.  This  extension  is  conceptually  straightforward  and  has  been 
proceeding  nicely.  A  three-dimensional  computational  mash  and  corres¬ 
ponding  storage  arrays  have  been  implemented.  Also,  numerous  convective, 
diffusive,  and  source  terms  involving  ©-derivatives  have  been  appropriately 
incorporated  into  each  of  the  governing  equations.  In  addition,  a  3-D 
iteration  scheme  has  bean  included. 

In  accordance  with  objectives  rsgarding  the  minimization  of  compu¬ 
tational  expense,  numerous  execution  cost  controls  have  been 
implemented.  As  previously  mentioned,  a  major  advancement  in  this 
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regard  is  the  incorporation  of  the  QUICK  differencing  scheme  in  an 
original*  and  stabilizing  manner*  which  has  dramatically  reduced  the 
cost  (via  reduced  grid  density)  required  for  necessary  solution  accuracy. 

This  scheme  was  found  to  yield  a  56*  execution  cost  reduction  for  a 
concentric-rotor  (2-D)  case  and  will  certainly  provide  considerably 
greater  reductions  for  the  all-important  eccentric-rotor  (3-D)  cases 
of  primary  interest  to  engine  manufacturers. 

i 

Several  other  coat  reducing  features  are  aimed  at  improving  the  con¬ 
vergence  rate  of  the  Iterative  solution  procedures.  Use  of  the  well-known 
Cyclic  Tri-Diagonal  Matrix  Algorith  as  well  as  the  SIMPLER  solution 
epproaoh  allow  fast  convergence.  Also*  an  advanced  data  Input/output 
sequence  has  already  provided  substantial  savings.  It  permits  writing 

to  and/or  reading  from  disk/tape  storage  the  currant  iteration  level 

\ 

values  at  any  desired  interval.  This  allows  the  user*  in  the  batch 
operating  mode,  to  monitor  solution  progress,  etc.*  and  adjust  the 
convergence  rate  parameters  and  then  resume  program  execution. 

F.  Summary 

Second  year  progress  in  developing  the  3-D  eccentric-rotor  computer 
code  included) 

(a)  incorporation  and  evaluation  of  the  QUICK  differencing 
schema  to  alleviate  false  diffusion  errors  and  reduce 
execution  costs, 

(b)  evaluation  of  the  2-D  concentric-rotor  version  via  comparison 
with  preliminary  measurements  of  compressible  flow  through 

an  abrupt  pipe  expansion. 
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(o)  brief  exploratory  prediction!  of  compressible  flow  within 
a  concentric-rotor  labyrinth  seal,  and 
(d)  partial  completion  of  extending  the  highly  successful 

concentric-rotor  code  to  the  required  eccentric-rotor  (3-D) 
version,  including  numerous  execution  cost  saving  features. 

G.  Graduate  Student  Development 

Mr.  Steve  Sobolik  continued  employment  during  the  second  year  of 
this  project  until  June  1,  1984.  He  received  a  B.S.  degree  (Cum  Laude) 
in  Mechanical  Engineering  from  TAMU  in  May  1982.  He  designed  and 
developed  software  for  Mitre  Corp.  (a  NASA  contractor)  during  three 
summers.  Steve  has  been  responsible  for  moat  of  the  computational 
development  work  on  this  project.  He  is  an  Innovative  individual  and 
produces  quality  work,  Ha  was  interviewed  by  numerous  AFRAPT  companies, 
and  received  an  offer  from  Garrett  Gas  Turbine  Engine  Co,  However,  ha 
accepted  a  more  lucrative  offer  elsewhere. 

Mr.  Steve  Hansel  began  employment  on  June  1,  1984,  having  received 
a  B.S.  degree  (Cum  Laude)  in  Nuclear  Engineering  from  TAMU  in  May  of  19S4. 
Mr.  Hemal  is  progressing  admirably  and  is  a  self-motivated,  Innovative 
individual., 

H.  Publications 

The  manuscript,  "Prediction  of  Incompressible  Flow  in  Labyrinth 
Seals,"  was  reviewed  by  the  ASME  Journal  of  Fluids  Engineering.  The 
manuscript  will  be  revised  slightly  to  comply  with  the  reviewers' 
suggestions. 
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I.  Conference  Presentations 

A  paper  consisting  of  ths  prsllainsry  prediction  test  csss  of  the 
sudden- expans ion  combustor  flow  field  entitled:  "Prediction  of  Subsonic 
Air  Flow  Through  a  Rocket/Ramjet  Combustor"  is  accepted  for  presentation 
at  the  AXAA  23rd  Aerospace  Sciences  Meeting.  It  will  be  held  January  14 
through  17,  1985  in  Reno.  Nevada.  The  authors  are  D.  L.  Rhode  and 
S.  R.  Sobollk.  Another  paper  consisting  of  the  brief  exploratory  pre¬ 
diction  of  compressible  flow  in  a  concentric-rotor  labyrinth  seal  is 
accepted  for  presentation  at  the  1985  ASMR  International  Gas  Turbine 
Conference.  The  title  is:  "Simulation  of  Subsonic  Flow  Through  a 
Generic  Labyrinth  Seal."  The  conference  is  being  held  March  18  through 
21,  1985  in  Houston,  Texas.  In  addition,  an  extended  abstract  for  a 
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SUMMARY 

i  i  i  1 

i  * 

i 

THa  baaio  Aquations  ara  dsrivad  for  oompresalble  flow  in  a 
labyrinth  aaal.  The  flow  is  assumed  to  be  completely  turbulent  in  the 
olroumferential  dlreotlon  where  the  frlotlon  factor  la  determined  by 
the  Blaaiua  relation.  Linearised  xeroth  and  first-order  perturbation 
equations  are  developed  for  small  motion  about  a  oentered  position  by 
an  expansion  in  the  eooentriolty  ratio.  The  soroth-order  pressure 
distribution  la  found  by  satiafylng  the  leakage  equation  while  the 
olroumferential  veloolty  distribution  is  determined  by  satisfying  the 
momentum  equation.  The  first  order  equations  are  solved  by  a 
separation  of  variable  solution.  Integration  of  the  resultant  pressure 
distribution  along  'and  around  the  seal  defines  the  reaotion  force 
developed  by  the  seal  and  the  corresponding  dynamio  coefficients.  The 
results  of  this  analysis  are  oompared  to  published  test  results. 

i 
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INTRODUCTION 

» 

The  problem  of  self  excited  vibration  In  turbonaohlnary  du«  to* 
labyrinth  aaala  haa  lad  to  tha  development  of  nany  analyaaa  whloh 
gttempt  to  nodal  tha  phyaloal  phanonanon  so  that  tha  problem  oan  ba 
better  understood  and  bharafora  aolvad.  Tha  shortcoming  with  tha 

analyaaa  whloh  hava  baan  praaantad  to  data  la  that  thay  ana  difficult 

* 

to  understand  and  require  limiting  aaaumptlona  auoh  aa  Ignoring  the 
araa  ‘  derivative  in  tha  olroumfarantial  diraotion,  aaauming  that  tha 

frlotion  factor  la  tha  aame  for  all  aurfaaea,  and  aaauming  that  tha 

* 

flow  ooaff Iolanta  la  oonatant  along  tha  seal.  Thaaa  aaaumptlona  may  ba 
of  aoma  uaa  mathematioally,  but  do  vary  little  for  tha  underatandlng  of 
tha  phyaloal  ooourrenoe.  Tha  firat  atapa  toward  analyala  of  thla 
problem  ware  taken  by  Alford  DL  who  neglootad  olroumfarantial  now 
and  Spur k  at  al.  [2]  who  naglaotad  rotation  of  tha  shaft.  Vanoa  and 
Murphy  [3]  extended  tha  Alford  analyala  by  introducing  a  more  raallstio 
aaaumption  of  ohokad  flow.  Xoatyuk  [4]  performed  tha  firat 
oomprehenalve  analyala,  but  failed  to  lnoluda  tha  change  in  araa  due  to 
eooentrioity  whloh  la  reaponsible  for  tha  relationship  between 
oroaa-Goupled  foroaa  and  parallel  rotor  diaplaoementa.  Xwataubo  C5,  6] 
refined  tha  Koatyuk  modal  by  lnoludlng  tha  time  dependency  of  araa 

t 

ohanga,  but  ha  naglaotad  tha  araa  darlvatlva  In  tha  olroumfarantial 
diraotion.  Kurohaahi  [7]  inoorporatad  dependency  of  tha  flow 
ooaff lolent  on  eooentrioity  Into  hie  analysis,  but  assumed  that  tha 
oiroumferential  veloolty  In  eaoh  oavity  was  tha  same. 

Tha  analysis  presented  hare  includes  tha  variation  of  tha  araa  In 
tha  olroumfarantial  diraotion  due  to  aooantrloity  and  incorporates  as 


“2- 


many  of  the  physical  phenomena  in  the  flow  field  as  was  thought 
neoaaaary  to  produo*  an  adequate  result.  The  main  purpose  of  this 
paper  is  to  present  a  unified  and  comprehensive  derivation  of  a  reduoed 
set  of  equations  and  a  new  solution  format  for  those  equations.  The 
results  of  this  analysis  are  oompared  the  published  test  results  of 
Waohter  and  Benokert  [8,  9,  10]. 

NOMENCLATURE 

Ai  ;  Cross  sootional  area  of  the  oavity  (L**2)i  defined  in  Eq.  (l.a) 

B*  Height  of  labyrinth  aeal  atrip  (L)i  defined  in  figure  (1) 
t  Dlreot  damping  coefficient  (Ft/L) 

Cr  Nominal  radial  olearanoe  (L)j  defined  in  figure  (1) 

Dh  Hydraullo  diameter  of  oavity  (L)i  introduced  in  Eq.  (3) 

H  Looal  radial  olearanoe  (L) 

K  Dlreot  stiffness  coefficient  (F/L) 

L  Pitoh  of  aeal  strips  (L)j  defined  in  figure  (1) 

NT  Number  of  seal  strips 

NC»NT"1  Number  of  oavltlea 
P  Pressure  (F/L*) 

R  Oaa  constant 

Ra  Radius  of,  seal  (L)j  defined  in  figure  (1) 

T  Temperature  (T) 

Ru  Surfaoe  veloolty  of  rotor  (L/t) 

Vi  Average  veloolty  of  flow  in  circumferential  direction  (L/t) 
a,b  Radial  seal  displacement  oomponents  due  to  elliptical  whirl  (L)j 
defined  in  Eq.  (13) 

ar  Dimension! eos  length  upon  which  shear  stress  aots  on  rotor 


as  Dimensionless  length  upon  which  shear  stress  aots  on  stator 


o  Croaa  ooupled  damping  ooeffloient  (Ft/L){  in  Eq,  (18) 
k  Croaa  oouplad  stiffness  coefficient  (F/L)j  in  Eq,  (18) 

A  Leakage  maaa  flow  rate  per  circumferential  length  (M/Lt) 

mr,  nr,  ms,  na  Coefficients  for  Blaaiua  relation  for  friction 
faotort  defined  in  Eq.  (3) 
t  Time  (t) 

w  Shaft  angular  velocity  (1/t) 

p  Density  of  fluid  (M/L»*3) 

v  Kinematio  viaooaity  (L**2/t) 

c  -  e/Cr  Eooentrioity  ratio 
Y  ■ '  Ratio  of  apeoif lo  heats 

Kq  Dlmansionlesa  crosa-ooupled  stiffness  parameter j  defined  in  Fqs.  (2 Y) 
Eq  Dimenaionleaa  entry  swirl  parameters  defined  in  Eqs.  (27) 

Subscripts 

o  Zeroth-order  component 

i  i-th  chamber  value 

1  First-order  component 

x  X-direotlon 

y  Y-dlreotlon 

r  Reservoir  value 

a  Sump  value 


M- 


PROCEDURE 


The  analysis  presented  here  is  developed  for  the  see-through  type 
of  labyrinth  seal  shown  in  figure  1,  The  continuity  and  momentum 
equations  will  be  derived  for  n  single  cavity  control  volume  as  shown 
in  figures  2,  3.  and  5.  A  leakage  model  will  be  employed  to  aooount 
for  the  axial  leakage.  The  governing  equations  will  be  linearized 
using  perturbation  analysis  for  small  motion  about  a  centered  position. 
The  seroth-order  continuity  and  momentum  equations  will  be  solvud  to 
determine  the  steady  state  pressure  and  velooity  for  eaoh  oavity.  The 
first-order  continuity  and  momentum  equations  will  be  reduoed  to 
linearly  independent,  algebraic  equations  by  assuming  an  elliptioal 
orbit  for  the  shaft  and  a  corresponding  harmonic  response  for  the 
pressure  and  velooity  perturbations.  The  foroe  and  force  ooefflolents 
for  the  seal  are  found  by  integration  of  the  first-order  pressure 
perturbation  along  and  around  the  shaft. 

ASSUMPTION 

1)  Fluid  is  oonsidered  to  be  an  ideal  gas. 

2)  Pressure  variations  within  a  ohamber  are  small  compared  to  the 
pressure  difference  aoross  a  seal  strip. 

3)  The  frequency  of  aoouatio  resonanoe  in  the  oavity  is  muoh  higher 
than  that  of  the  rotor  speed. 

*0  Added  mass  terms  are  negleoted. 

5)  The  eocentrioity  of  the  rotor  is  small  compared  to  the  radial  seal 
clearance,. 
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6)  In  the  determination  of  the  shear  stresses  in  the  oiroumferential 

direction,  the  axial  component  of  velooity  is  neglected. 

7)  The  contribution  of  shear  stress  to  the  stiffness  and  damping 

ooeffioienta  is  negleoted. 

GOVERNING  EQUATIONS 


Continuity  Equation 

Referring  to  the  control  volume  in  figures  2  and  3,  the  continuity 

equation  for  the  control  volume  shown  is: 

*  » 

3  3  PiViAi 

(piAi)  +  *  +  "  0  (1) 

where  the  transverse  surfaoe  area  A*  is  defined  byj 

Ai  -  <Bi  +  ♦  Bi+1  +  H1+1)  L1/2  Cl .a) 

Momentum  Equation 

4 

The  momentum  equation  (2)  is  derived  using  figures  4  and  5  which 
show  the  pressure  foroes  and  shear  stresses  acting  on  the  oontrol 

volume.  This  equation  includes  the  area  derivative  in  the 
oiroumferential  direction,  which  was  negleoted  by  Iwatsubo  [3,  6]. 

dpViAi  2pViAi  3Vi  3Vai  Hi  V^i  3p 

3b  -  Rs  30  Rs  30  Rs  30  (2) 

Ai  3Pi 

-  mi+1Vi  -  miVi-!  -  "  “  “  +  TriariLi  "  TslaolLl 


where  ar  and  as  are  the  dimensionless  length  upon  which  the  shear 
stresses  aot  and  are  defined  for  teeth  on  rotor  by 


and  for  teeth  on  the  stator  by 

as*  -  (2B*  Li)/Li  arj  -  1 

Blaaiua  [11]  determined  that  the  shear  stresses  for  turbulent  flow 

In  a  smooth  pipe  could  be  written  as 

t  mo 

S  •  i.  pUffl  no  ^UniDhj 

where  Urn  is  the  mejan  flow  velocity  relative  to  the  surface  upon  whioh 
the  shear  stress  is  aoting.  The  constants  mo  and  no  oan  be  empirically 
determined  for  a  given  surfaoe  from  pressure  flow  experiments. 
However,  for  smooth  surfaces  the  ooeffioiants  given  by  Yamada  [12]  for 
turbulent  flow  between  annular  surfaoea  arei 

i 

mo  ■  -0.25  no  ■  0.079  • 

Applying  Blaslus'  equation  to  the  labyrinth  rotor  and  stator  surfaces 
yields  the  following  definitions  for  the  rotor  and  stator  shear 
stresses. 


where  Dh^  is  the  hydraulic  diameter  defined  by 


Separate  parameters  (ms,  ns),  (mr,  nr)  for  the  stator  and  rotor, 
respectively ,  aooount  for  different  rotor  and  stator  roughnesses. 


If  Eq.  tl)  times  the  circumferential  velooity  is  now  subtracted 
from  Eq.  (2),  the  following  reduotd  form  of  the  momentum  equation  is 
obtained! 


1  +  Sr*  <*1  ’  V "  u‘ " 1 


-A, 


1  9P 


T,l*,iLi  <5) 


+  -Wfr. 


In  order  to  reduoe  the  numbers  of  variables,  all  of  the  density  terms 
are  replaoed  with  pressure  terms  using  the  ideal  gas  law. 

Pi  “  pRT  (6) 

» 

Furthermore,  in  order  to  make  the  perturbation  analysis  easier,  the 
following  substitution  is  made  in  the  continuity  equation! 

.  .  *1+1-  *i 

*1+1  “  Ai  2ifi  • 

o 

Leakage  Equation 

To  aooount  for  the  leakage  mass  flow  rate  in  the  continuity  and 
momentum  equations,  the  leakage  model  of  Neumann  [13]  was  chosen.  This 
model  prediots  leakage  and  pressures  fairly  aoourately  and  has  a  term 
to  aooount  for  kinetio  energy  carryover.  However,  the  empirical  flow 
coefficient  relations  given  by  Neumann  were  discarded  in  favor  of  the 
equations  of  Chaplygin 


cm]  for  flow  through  an  orlfloe.  This  was  done  to  produoe  a  different 

V 

flow  ooeffiolent  for  suooeeding  oontraotions  along  the  seal,  aa  has  been 

*  • 

shown  to  be  the  oaae  by  Ggll  [153.  The  form  of  the  model  1st 


*1  "  yliy2  V  ^  *1-1  "  ?1 


KT 


(7) 


where  the  kinetic  energy  oarryover  ooeffiolent  pa  is  defined  for 
straight  through  seals  ast 


where 

'J  ■  1  -  (1  ♦  16.6  Cr/L)-2 


NT 

ThWTT 


and  is  unity,  by  definition,  for  interlooking  and  oombination  groove 
seals.  The  flow  ooeffiolent  is  defined  ast 


it 

ir+2-50.j+2a£ 


where  a^ 


-  1 


For  ohoked  flow,  Fliegner's  formula  [16]  will  be  used  for  the  last  seal 
strip.  It  is  of  the  formt 


PERTURBATION  ANALYSIS 


For  oavity  1,  the  continuity  Eq.  (1),  momentum  Eq.  (5)  and  leakage  Eq. 
(7)  are  the  governing  equations  for  the  variables  V*i  P**  A**  A 
perturbation  analysis  of  these  equations  is  to  be  developed  with  the 
eooentrloity,  ratio,  c  -  e/Cr,  seleoted  to  be  the  perturbation 

,  i 

parameter.  The  governing  equations  are  expanded  in  the  perturbation 

t 

variables 

Pi  “  P0i  ♦  t  Pii  Hi  -  Crl  ♦  c  Hi 

Vi  *  Voi  +  c  VU  Ai  -  A0  ♦  e  KHi 

where  c  -  e/Cr  is  the  eooentrloity  ratio.  The  zeroth-order  equations 
define  the  leakage  mass  flow  rata  and  the  olroumfarentlal  velootty 
distribution  for  a  oentered  position.  The  first-order  equations  define 
the  perturbations  in  pressure  and  olrounferential  velooity  due  to  a 
radial  position  perturbation  of  the  rotor.  Strlotly  speaking,  the 
results  are  only  valid  for  small  motion  about  a  oentered  position. 

Zeroth-Order  Solution 

The  zeroth-order  leakage  equation  is 

“i  +  ,1  "  *i  “  n»o  (9) 

and  is  used  to  determine  both  the  leakage-rate  m0  and  pressure 
distribution  for  a  oentered  position.  The  leakage  rate  is  determined 
using  either' Eq.  (7)  or  Eq.  (6),  depending  on  the  operating  conditions. 
To  determine  if  the  flow  is  ohoked  or  not,  assume  that  the  pressure  in 


the  last  oavity  is  equal  to  the  orltioal  pressure  for  ohoklng.  Using 
this  pressure,  find  the  leakage  from  Eq.  ‘(8)  and  then  use  Eq.  (7)  to 
determine  the  reservoir  pressure  neoessary  to  produoe  this  oondition. 
Based  on  this  pressure,  a  determination  oan  be  made  whether  the  flow  is 
ohoked  or  not.  The  asaooiated  pressure  distribution  is  determined  by 

t 

employing  the  oaloulated  leakage,  along  with  a  known  boundary  pressure, 
and  solving  Eq.  (7)  sequentially  for  eaoh  oavity  pressure  F0l. 

The  seroth“>order  oiroumferentlal-momentum  equation  is 
m0  (V0i  -  VoM)  ■  (trio  ari  "  ?slo  i8i)  l1  "  1*  2,  • .  .NC , 

Prom  oaloulated  pressures,  the  densities  oan  be  oaloulated  at  eaoh 
oavity  from  Eq.  (6),  and  the  only  unknowns  remaining  in  Eq.  (10)  are 

the  oiroumferential  velocities  V0i*  Given  an  inlet  tangential 

«  > 

velooity,  a  Hewton-root-findlng  approaoh  oan  be  used  to  solve  Eq.  (10) 
for, the  1-th  velooity,  one  oavity  at  a  time.  This  is  done  starting  at 

1  i 

the  first  oavity  and  working  down  stream. 


First-Order  Solution 

The  governing  first-order  equations  (11,  12),  define  the  pressure 

i 

and  velooity  fluotuations  resulting  from  the  seal  olearanoe  function. 
The  oontinuity  equation  (11)  and  momentum  equation  (12)  follow: 

0uTr+0li'%,Tr+l!iiT?  -TT+°31  fli  +  04l  ru-l  +  GJ1 

•  8Hi  „ 

"  "  °6i  Hi  **  °2i  “Tt  “  G2i  Rs  36 


r  >Vli  XUVol  !!u  +  ^oi^U  +  x  y 
XllT T  +  T5—  ~+  Rb  36  '  2iVli 


Ao  Vli-1  +  X3i  P11 


*  *41  *11-1  -  X5i  H1 


where  the  Xj's  and  G^'a  are  defined  in  Appendix  A.  If  the  shaft  center 
moves  In  on  elllptloal  orbit,  then  the  aeal  ol ear anoe  function  oan  be 
defined  ast 


I 

1 

I 


1 

E§ 

1 

1 

l« 


cH|  ■  -a  ooswt  ooefi  -b  si nut  sln6  (13) 

■  -a  [oos  (6-uit)  +  oos  (e+wt)3  -  b  Coos  (8-wt)  *  oos(8+wt)3 

1  2 

The  pressure  end  velocity  fluctuations  oan  now  be  stated  in  the 
associated  solution  format) 

p1i  -  Poi°os(8+ut)  +  PJiSin(8+«t)  ♦  PoioosCe-mt)  ♦  PgiSlnte-ut)  (14) 

i 

VU  -  vSj,ooa(e+wt)  +  VjiSin(8+ut)  +  V5ieos(8-wt)  +  V^sinO-mt)  (15) 

Substituting  Eqs.  (13),  (14),  and  (15)  into  Eqs.  (11)  and  (12)  and 

grouping  like  terms  of  sines  and  oosines  (as  shown  in  Appendix  B) 

( 

eliminates  the  time  and  theta  dependency  and  yields  eight  linear 
algebraic  equations  per  oavity.  The  resulting  system  of  equations  of 
the  i-th  oavity  is  of  the  form) 

CAj-i ]  (Xi-i)  +  [Aj]  (Xj,)  +  [Aj+13  (Xj+1 )  -  »  (B^)  +  b  (ci)  06) 

c.  7. 

where 

(Xi-i)  -  (Pji-i,  Pol-1 •  psi-1i  Poi-1*  vsi-1»  vJi-,.  vJi-i,  V0l-i)T 

«!>  -  (p;t,  r*oi,  pii,  p^i, v*!.  »;t,  v;i)T 

Mi  ^  i^a  m 

(Xi+1 )  -  (P8i+1*  Pcl+1  •  psi+1.  poi+1*  vsi+1»  vci+1*  vci+1*  vsi+1 ) 

The  A  matrio.es  and  column  veotors  B  and  C  are  given  in  Appendix  B.  To 
use  Eq,  (16)  for  the  entire  seal  solution,  a  system  matrix  must  be 
formed  whioh  is  block  tridiagonal  in  the  A  matrices,  The  size  of  this 
resultant  matrix  is  (cNC  X  cNC)  sinoe  pressure  and  velooity 
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perturbations  at  the  inlet  and  the  exit  are  assumed  to  be  zero*  This 
system  is  easily  solved  by  various  linear  equation  algorithms,  and 
yields  a  solution  of  the  forms 


PSi  -  a  Fast  ♦  b  F6si 

t  * 

psi  ■  a  Fasi  +  b  Fbsi 

c  c  (17) 

Pol  “  a  Faoi  b  Ffcoi 

i  c 

pol  ■  a  P'ioi  *  b  Fboi 

.7  e 

DETERMINATION  OF  DYNAMIC  COEFFICIENT 


The  foroe-motlon  equations  for  a  labyrinth  seal  are  assumed  to  be 
of  the  forms 


l 

Fx 

i  r*  ki  | 

X)  [*C  < 

>1  lxl 

I 

"  1 

M 

<18> 

1 

py 

1  L-k  kJ  1 

Y '  L-o  ( 

)J  lYl 
• 

The 

solution  of  Eq. 

(18)  for  the 

stiffness 

and  damping  ooeffioients  is 

the 

objeotive  of 

the  ourrent 

analysis . 

For  the  assumed  elllptioal 

orbit  of  Eq.  (13).  the  X  and  Y  components  of  displaoement  and  velocity 
are  defined  ass 

X  ■  a  ooawt  X  -  -am  sinust 

i 

Y  ■  b  sinust  Y  ■  bus  oosust 
Substituting  these  relations  into  (18)  ylsldst 

+FX  "  -Ka  ooaust  -  kb  si  nut  +  Caus  sinust  -  cbui  oosust  (19) 

•s 

+Fy  *  -ka  ooaust  -  Kb  sinust  -  oaus  sinust  -  Cbus  oosust 


I 

I 


I 

I 

I 

B 


|| 


N 

19 


I 


Redefining  the  forces,  Fx  and  Fy,  aa  the  following! 

Fx  *  Fxo  0Q3<ut  +  Fxa  flintut  (20) 

Fy  m  FyQ  ooaut  ♦  Fyg  slnwt 

and  aubatltutlng  baok  Into  (19)  yields  the  following  relational 

-Fxo  •  Ka  ♦  obw  -Fxa  -  -CaX  ♦  kb  (21 ) 

-Fyo  «  ka  ♦  Cbu  ~Fys  ■  Kb  ♦  caw 

i 

The  X  and  X  oomponenta  of  foroe  can  be  found  by  Integrating  the 

i  • 

preaaure  around  the  seal  aa  follows: 


NC 

2v 

1 

Fx  •  -Rac 

I 

/ 

Pli  Li  oosO  dfl 

(22) 

i-1 

0 

NC 

2v 

Fy  ■  Rats 

I 

/ 

Pu  Li  alnB  dB 

(2B) 

i-1 

0 

Only  one  of  these  oomponenta  needs  to  be  oxpandod  in  order  to  determine 
the  dynamic  ooeffioients.  For  this  analysis,  the  X  oomponent  waa 
ohoaen,  Substituting  Gq.  (I1))  into  (22)  and  integrating  yieldat 

Fx  ■  -ctRa  i  Lj  C(Pji  -  pJi)  einwt  ♦  (P0i  ♦  Poi)  ooewt]  (2i|) 
1-1 

substituting  from  Eq.  (17)  and  (19)  into  Eq,  (2*0  and  equating 
ooeffioients  of  si  nut  and  ooaut  yields! 


NC  .  +  -  +  - 

Fxa  *  -irRo  j  Ca(Faaj  ~  Faai )  +  b  (Fpai  ~  Pbai^-9 

i-1 

»»  ■  "S  W  Ca(Fi„i  +  f;o1)  *  I>  <r£0l  *  FboiU 

■,  1-1 


(25) 


Equating  the  definitions  for  Fxa  and  Fxo  provided  by  Eqs.  (21)  and  (25) 
and  grouping  like  terns  of  the  linearly  Independent  ooeffloients  a  and 
b  yields  the  final  solutions  to  the  stiffness  and  danplng  coefficients! 


k  ■  »R  l  (Fbai  "  pbsi  ^1 
1-1 


-wBa  NC  +  1  - 

C  -  w  l  (f*si  “  Fasi' 

i-1 

wRa  NC  +  - 

c  ■  u>  (Fbol  +  Pbei  L1 

SOLUTION  PROCEDURE  SUhMARY 

In  review,  the  solution  procedure  uses  the  following  sequential  stepsi 

a)  Leakage  is  deternlned  from  Eq,  (7)  or  (8). 

b)  Pressure  distribution  is  found  using  Eq.  (7). 

o)  Veloolty  distribution  is  deternlned  using  Eq.  (10). 

d)  A  system  equation  is  formed  and  solved  using  the  oavlty 

Eq.  (16). 

e)  Results  of  this  solution,  as  defined  in  Eqs.  (17) .  are  Inserted  Into 

Eq.'  (26). 

RESULTS 

% 

To  oompare  the  present  analytic  solution  with  the  experimental 
results  of  Waohter  and  Beokert  [8,  9,  10],  the  following  dimensionless 
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parameters  are  introduced.  The  dimensionless  oroas- coupled  stiffness 
and  entry-swirl  parameters  are  defined  by  Waohter  and  Benckert  asi 


#  Cr  Kxy 
K(j  «  RsLW(Pr-Pa) 


e 


- 


0.5d»V^ 

(Pr-Ps)  I  0?5po^x 


(27) 


Ail  of  ths  results  presented  for  comparison  In  this  paper  ere  for  s 
seal  with  teeth  on  the  stator,  with  entry  swirl,  and  no  shaft  rotation. 
Although  Waohter  and  Benkert  published  results  for  shaft  rotation,  the 
data  for  the  operating  conditions  and  seal  geometry  were  insufficient 
for  use  In  this  study.  The  results  in  figures  7*  8,  and  9  are  from  [8] 
and  show  the  relationship  between  cross  coupled  stiffness  and  ths  entry 
.swirl,  for  a  seal  with  strips  on  the  stator  and  the  geometry  shown  in 
figure  6.  The  lino  shown  la  the  experimental  result  and  the  symbols 
are  the  results  from  this  analytical  modal.  These  figures  show  that 
the  modol  oompares  favorably  to  the  experimental  results  in  magnitude 
and  the  overall  trend  for  various  operating  oondltions,  Ths  figures 
also  show  that  the  model  does  not  yield  a  consistently  high  or  low 
result.  Instead,  the  model  tends  to  over  prediot  the  value  of  the 
stiffness  for  a  large  number  of  strips  and  under  prediot  stiffness  for 
a  small  number  of  seals.  This  trsnd  is  probably  dus  to  srrors  in 

i 

oaloulating  the  zeroth-order  pressure  distribution  using  the  leakage 
model . 

t 

The  results  in  table  (1)  are  from  [9,  10]  for  s  seal  with  strips 

on  the  stator.  The  results  show  the  effect,  of  change  in  seal 
parameters  suoh  ae  pitch,  number  of  teeth,  radius,  strip  height,  and 
olearanoe  on  the  oross  ooupled  stiffness.  The  model  aoouratoly  shows 
the  increase  in  cross-coupled  stiffness  due  to  decrease  in  olearanoe 
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and  deorease  in  atrip  height,  but  it  fails  to  remain  oonatant  Tor 

ehanga  of  pitoh  and  consistently  over  ea'imatea  the  oross-ooupled 

« 

•tiffneaa  for  the  larger  radlua  oaaea  by  about  26}. 

,  CONCLUSION 

A  olear  and  underatandable  analysis  utilising  reduced  equations 

has  .  been  presented  for  the  1  problem  of  oaloulatlng  rotordynamic 

ooeff loients  for  labyrinth  seals.  This  paper  was  developed  to  provide 

% 

a  leas  reatrlotlve  analysis  and  a  better  explanation  of  the  current 
analyses.  The  model  developed  gives  results  that  are  within  25%  of  the 
experimental  results  which  are  available.  However,  this  error  must  be 
balanoed  against  the  known  uncertainties  in  the  experimental  data. 
This  is  especially  Important  since  all  of  the  data  used  are  for  a 
nonrotating  shaft  and  the  only  lnflusnoe  on  the  oross  coupled  stiffness 
was  the  entry  swirl.  Although  Waohter  and  Bsnokert  published  data  for  a 
rotating  shaft,  the  data  were  not  sufficient  to  calculate  a  result. 
Also,  the  only  data  available  for  see-through  labyrinths  is  for  the 
type  with  strips  on  the  stator.  For  a  more  rigorous  test  of  this  and 
other  models,  more  complete  data  are  required  ovsr  a  wider  range  of 
parameters  for  different  seal  geometries.  Finally,  this  analysis  is 
only  oonsldered  valid  for  the  see-though  type  of  labyrinth  seal  since 
the  model  fared  very  poorly  in  comparison  with  experimental  results  for 
Interlooking  and  grooved  seal  data. 


-17- 


/IRA 


■»MV 


REFERENCES 


1.  Alford,  J.  S.,  "Protacting  Turbomaohinery  from  Self-Exoited  Rotor 
Whirl,"  Tranaaotlona  ASME  J.  of  Engineering  for  Power.  Ootober  1965, 
pp.  333-3^7 

2.  Spurk,  J.  H,,  and  Kaiper,  R.,  "Selbaterregte  Sohwingungen  b«i 
Turboraaaohinen  infolga  der  Labyrlnthatromung,"  Ingerleur-Arohive  M3, 

1974,  pp.  127-135. 

3.  Vanoa,  J.  M.,  and  Murphy,  B.  T.,  "Labyrinth  Sail  Effaota  on  Rotor 
Whirl  Stability,"  Inat.  of  Maohanioal  Engineer,  1980,  pp.  369-373. 

4.  Koatyuk,  A.  0.,  "A  Thaoratloal  Analysis  of  the  Aerodynamic  Foroea 
in  the  Labyrinth  Glands  of  Turboraaohinoa,"  Teploenergetloa,  19  (11)0,  , 
1972,  pp.  39-44. 

5.  Xwatsubo,  T.,  "Evaluation  of  Instability  Foroea  of  Labyrinth  Seala 
in  Turbines  or  Compressors,"  NASA  CP  2133  Proceedings  of  a  workshop  et 
Texas  A AM  University  12-14  May  1980,  Entitled  Rotordynamio  Instability 
Problems  in  High  Performance  Turbomaohinery,  pp,  139-167. 

6.  Xwatsubo,  T.,  Matooka,  N.,  and  Kawai,  R.,  "Flow  Induced  Force  and 
Flow  Pattern  of  Labyrinth  Seal,"  NASA  CP  2250  Proceedings  of  a  workshop 
at  Texas  AAM  University  10-12  May  1982,  Entitled  Rotordynamio 
Instability  Problems  in  High  Performance  Turbomaohinery,  pp.  205-222. 

7.  Kurohaahi,  M.,  Inoue,  Y« ,  Abe,  T,,  and  Fujikawa,  T.,  "Spring  and 
Damping  Coefficients  of  the  Labyrinth  Seal,"  Paper  No.  C2B3/80  delivered 
at  the  Seoond  International  Conferenoe  on  Vibrations  in  Rotating 
Maohinery,  The  Inst,  of  Meoh.  Engineering. 

8.  Waohter,  J.,  and  Benokert,  H.(  "Querkrafte  aua  Spaltdiohtungen-Elne 
mogliohe  Uroaohe  fur  die  Laufunruhe  von  Turbomaaohinen,"  Atomkernenergle 
Bd.  32.  1978,  Lfg.  4,  pp.  239-246. 

9.  Waohter,  J.,  and  Benokert,  H.,  "Flow  Induoed  Spring  Coefficients  of 
Labyrinth  Seala  for  Applications  In  Rotordynamio,"  NASA  CP  2133 
Proceedings  of  a  workshop  held  at  Texas  AAM  University  12-14  May  1980, 
Entitled  Rotordynamio  Instability  Problems  of  High  Performance 
Turbomaohinery,  pp,  189-212. 

10.  Benokert,  H,,  "Stromungsbedinte  Foderkennwerte  in 
Labyrlnthdichtungen,"  Doctoral  dissertation  at  Univeraity  of  Stuttgart, 
1980. 


11.  Blaslus,  H,,  "Forsohungoarb,"  Ing,-Wes,(  No  131,  1913. 

12.  Yamada,  Y,,  Trans.  Japan  Soc.  Meohanioal  Engineers,  Vol.  27,  No. 

180,  1961,  pp.  1257:  "  ' 


-18- 


13*  Neumann.  K.,  "Zur  Frage  der  Verwendung  von  Durohbliokdichtungen  im 
Dampfturblnenbau,"  Maaohlnenteohnlk ,  Vol.  1 3 »  196M ,  No  *t. 

1H.  Gurevieh,  M.  I..  The  Theory  of  Jets  In  An  Ideal  Fluid.  Pergamon 
Prase,  1966,  pp,  319-323. 

15.  Egli,  A.,  "The  Leakage  of  Steam  Through  Labyrinth  Glands,"  Trans. 
ASME,  Vol.  57,  1935,  pp.  115-122. 

16.  John,  E,  A,  James,  Oas  Dynamloa,  Wylie,  1979. 


m 

1 

m 


l 


r 

i  ->• 


i 

* 


m 

| 

* 

! 

I 


I 


XV  \ 
w 


CO 


Figure  4 
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Appendix  A 


Definition  of  the  First-order  Continuity 
and  Momentum  Equation  Coefficients 
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A1_1  MATRIX 


*1,2  "  a2,l  "  *3,4  “  *4,3  "  G4 

*5,2  "  *6,1  "  *7,4  "  *8,3  "  X4 

*5,6  "  *6,5  "  *7,8  "  *8,7“Ao 
The  remaining  element!  are  zero 


Appendix  B 


Separation  of  the  Continuity  and  Momentum  Equations 
and  Definition  of  the  System  Matrix  Elements 
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A  Comparison  of  Experimental  and  Theoretical  Results  for  Leakage, 
Pressure  Distribution,  and  Rotordynamic  Coefficients  for  Annular  Qaa 

Seals.  (Deoember  1984) 

Colby  Oran  Nicks,  B.S.,  Virginia  Polyteohnio  Institute  and  State 

University 

Chairman  of  Advisory  Committee:  Dr.  Dara  Childs 

This  thesis  conoerns  a  study  of  annular  gas  seals  which  is 
currently  in  progress  at  Texas  A&M  University.  A  brief  discussion  of 
the  importance  of  seal  behavior  in  rotordynamios  is  presented,  as  is  a 
review  of  current  annular  seal  theory.  An  outline  of  Nelson's 
analytical-computational  method  for  determining  rotordynamic 
coefficients  for  this  type  of  oom pr ess i ble- flow  seal  Is  included. 
Various  means  for  the  experimental  identification  of  the  dynamic 
coefficients  are  outlined,  and  the  method  employed  at  the  TAMU  teat 
facility  is  explained.  The  TAMU  test  apparatus  is  described,  and  the 
test  procedures  are  disoussed.  Experimental  results,  including  leakage, 
entranoe-loss  coefficients,  pressure  distributions,  and  rotordynamic 
ooeffioients  for  a  smooth  and  a  honeycomb  oonstant-olearanoe  seal  are 
presented  and  oompared  to  theoretical  results  from  Nelson's  analysis. 
The  results  for  both  seals  show  little  sensitivity  to  the  running  speed 
over  the  test  range.  Agreement  between  test  results  and  theory  for 
leakage  through  the  seal  is  satisfactory.  Test  results  for  direct 
stiffness  show  a  greater  sensitivity  to  fluid  prerotation  than 
predioted.  Test  results  show  that  the  deliberately-roughened  surface  of 
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the  honeycomb  seal  provides  improved  stability  versus  the  smooth  seal. 
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NOMENCLATURE 


■  Fourier  coef'f icienta  for  rotor  motion 


ns  ,ms 


nr  „mr 


direct  and  orosa-coupled  stiffness  ooeffloients  (FT/L) 


displacement  of  seal  rotor  from  centered  position  (L) 


direct  and  orosa-coupled  stiffness  coefficients  (F/L) 


entrance-loss  coefficient 


direot  and  orosa-coupled  added-maas  coefficients  (M) 


fluid  mass  flow  rate  (M/T) 


stator  Hire  coefficients 


rotor  Hirs  coefficients 


fluid  pressure  (F/L*) 


seal  radius  (L) 


2pUC/u  ■  nominal  axial  Reynolds  number 


mean  fluid  flow  velocity  (L/T) 


radial  seal  displacements  (L) 


ratio  of  specific  heats  for  air 


«o  >  Cr  -  equilibrium  eocentrioity  ratio 


fluid  density  (M/L‘) 


Fanning  friocion-factor 


fluid  shear  stress  (F/La) 


shaft  angular  velocity  (1/T) 


shaft  preoessional  velocity  (1/T) 


fluid  viscosity  (FT/L* ) 


INTRODUCTION 


With  turbomaohinery  design  trends  tending  toward  increased  speeds 
and  loadings,  lighter  weight,  and  reduced  olearanoes  between  rotating 
and  stationary  parts,  oonsiderable  ooncern  with  instability  and 
synchronous  response  has  arisen.  Synchronous  response  refers  to 
vibration  of  the  turbonaohine  rotor  assembly  at  a  frequency  coincident 
with  the  rotational  speed.  Characteristically,  the  vibration  amplitude 
inoreas.es  to  a  maximum  at  eaoh  orltioal  speed  (coincidence  of  the 
running  speed  with  a  rotor's  damped  natural  frequency),  and  then 
decreases  to  a  relatively  steady  level.  Operation  of  turbomaohines  at 
rotational  speeds  above  any  of  the  critioal  speeds  requires  the  rotor  to 
traverse  them  during  start-up  and  shut-down.  Therefore,  in  order  to 
limit  the  peak  synchronous  vibration  levels,  the  maohine  designer 
aspires  to  jntroduoe  damping  into  the  rotor  system. 

In  contrast  to  synchronous  vibration,  "unstable"  or  "self-excited" 
motion  is  typloally  subaynohronous .  This  motion  takes  the  form  of 
whirling  of  the  rotor  shaft  at  a  natural  frequency  less  than  the 
rotational  speed.  The  exoiting  force  for  this  whirling  motion  is  a 
tangential  force  aoting  on  the  rotor  due  to  some  fluid  or  friction 
meohanJ.sm.  This  vibration  often  ocours  with  large  amplitudes  which 
sustain  or  grow  as  running  speed  increases.  At  beat,  this  self-exoited 
whirling  prevents  further  speed  increases)  at  worst,  it  results  in 
damage  to  or  oatastrophio  failure  of  the  equipment.  One  of  the 
rotordynamio  force  mechanisms  which  plays  a  role  in  self-exoite^  rotor 
motion  and  synchronous  response  is  that  of  the  foroes  developed  by 
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annular  seals.  Until  recently,  meat  investigations  of  annular  seals  in 
turbomachinery  have  been  concerned  wi-th  reduoing  the  leakage  of  the 
working  fluid  through  the  seal  (i.e.,  improving  the  sealing  effect). 
However,  reoent  experiences  have  shown  that  forces  developed  by  these 
seals  oan  have  considerable  lnfluenoe  on  the  stability  and  synchronous 
response  of  rotating  machinery.  Black  et  al.  [1-33  have  demonstrated 
the  oritioal  effects  that  forces  developed  by  neck-ring  and  interstage 
seals  have  on  the  rotordynamic  behavior  of  pumps.  Also,  stability 
difficulties  with  the  Space  Shuttle^  Main  Engine  (SSME)  high-pressure 
fuel  turbopump  [M]  have  prompted  further  research  into  these  forces 
developed  by  liquid  seals. 

Experiences  have  shown  that  various  gas  seal  configurations  oan 
have  similar  influences  on  the  rotordynamic  behavior  of  turbomachinery . 
In  the  high-pressure  oxygen  turbopump  of  the  SSME,  for  example,  initial 
vibration  problems  were  remedied  by  changing  the  turbine  interstage  seal 
from  a  atepped-labyrinth  configuration  to  a  convergent  taper  seal  with  a 
honeycomb  stator  and  a  smooth  rotor[5].  A  lack  of  experimental  data  to 
completely  explain  this  and  other  gas  seal  behavior  makes  obvious  the 
need  for  research  in  this  area. 

The  purpose  of  this  report  is  twofold.  It  desorlbes  the  test 
facility  and  initial  test  program  developed  to  experimentally  measure 
the  fluid  forces  induced  by  annular  gas  seals,  and  it  provides  a 
comparison  of  theoretically  predicted  and  experimentally  obtained  data 
for  smooth  and  honeycomb  soals.  The  leakage  of  the  working  fluid 
through  the  seal,  the  pressure  gradient  along  the  seal  length,  entranoe 
pressure-loss  data,  and  rotordynamic  coefficients  provide  a  basis  for 
comparison.  A  short  discussion  on  seal  theory  is  included,  and  various 


ANNULAR  SEAL  ANALYSIS  REVIEW 


Aa  related  to  rotordynamioa,  seal  analysis  has  the  objective  of 
determining  the  reaotion  forces  acting  on  the  rotor  arising  from  shaft 
motion  within  the  seal.  Due  to  similarities  between  plain  Journal 
bearings  and  annular  seals,  seal  analysis  is  generally  based  on 
governing  equations  whioh  have  previously  been  developed  for  bearings. 

Annular  seals  and  plain  bearings  are  geometrically  similar,  but 
seals  typically  have  radial  olearanoe-to-radius  ratios  on  the  order  of 
0.005,  versus  Cr/R  ratios  of  0.001  for  bearings.  Due  to  seal  clearances 
and  pressure  differentials,  fully-developed  turbulent  flow  normally 
exists.  Also,  seals  are  nominally  designed  to  operate  in  a  centered 
position.  Journal  bearings,  on  the  other  hand,  have  operating 
eaoentriolties  whioh  vary  with  running  speed  and  load.  Therefore,  moat 
of  the  rotordynamio  work  for  bearings  has  been  done  to  determine  dynamic 
coefficient  versus  eccentricity  relationships. 

Two  linearized  seal  models,  expressed  in  terms  of  dynamic 
coefficients,  have  been  suggested  for  the  motion/reaotion-foroe 
relationship.  For  small  motions  of  the  rotor  about  an  arbitrary 
position  In  the  seal,  as  shown  in  Fig.1,  the  relation  oan  be  written 

'  i 

FX  KXX<« o>  KXY<« o>  x  "CXX<c0>  CXY<£o>  X  Mxx<co>  MXY<£o>  ft 

:W  +  ’+  1  (  1  ) 

FY  KYX<*o)  KYY^co^  y  CYX^co^  CYy(®o)  ft  MYX^eo5  MYy(£o)  ft 

where  the  dynamic  coefficients  (KXX,  KYY ,  CXX,  CYY ,  MXX,  MYY)  and  ( KXY , 
KYX,  CXY,  CYX,  MXY,  MYX)  represent  the  "direct"  and  "cross-coupled" 
stiffness,  damping,  and  added-mass  terms,  respectively .  These 


Fig.  1  Small  motion  of  a  aaal  rotor  about  an  accantric  position; 
u»  la  tha  rotor  spin  spaed, A  ia  tha  processional  orbit 
frequency. 


Fig.  2  Small  motion  of  a  aaal  rotor  about  a  cantered  position; 
w  ia  tha  rotor  spin  speed, A  is  tha  precassional  orbit 
frequency. 
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coefficients  are  functions  of  the  equilibrium  eocentrioity  ratio 
c0  -  a0  /  Cp,  where  the  eooentrioity  ratio  c0  equals  the  displacement 
(a0)  of  the  rotor  from  the  centered  position  divided  by  the  nominal 
radial  clearance  (Cr).  The  term  "cross-coupled"  refers  to  the  coupling 
effect  exhibited  by  the  off-diagonal  terms;  specifically,  motion  in  one 
plane  introduces  reaotlon  forces  in  an  orthogonal  one.  These 
cross-coupled  terms  arise  from  the  fluid's  oiroumferential  velooity 
component,  and  show  a  strong  dependency  on  both  the  magnitude  and 
direction  (with  respect  to  rotor  rotation)  of  the  velooity.  This 
oiroumferential  velocity  component  may  arise  from  the  prerotation  of  the 
fluid  as  it  enters  the  seal  due  to  some  rotating  element  upstream,  or  it 
may  develop  as  the  fluid  passes  through  the  seal,  with  rotor  shear 
foroes  "dragging"  the  viscous  fluid  around  its  periphery.  The 
oross-couplad  stiffness  term  usually  produces  a  destabilizing  foroe 
component,  and  therefore  is  of  considerable  interest.  The  oross-coupled 
damping  and  added-mass  terms  are  generally  much  leas  influential  than 
the  oroaa-ooupled  stiffness  term  with  respect  to  stability.  For  no 
fluid  rotation,  these  oroas-ooupled  terms  are  zero. 

The  second  linearized  seal  model  applies  for  small  motions  of  the 
rotor  about  a  oentered  position  in  the  seal,  as  shown  in  Fig.  2  .  This 
model  oan  be  expressed 


K 


(2) 


where  the  dynamic  coefficient  matrioes  are  skew-symmetric. 

Theoretical  work  on  annular  seals  has  been  done  for  both 
incompressible  and  compressible  fluids.  Black  et  al .  [6]  have  developed 
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analytical  "ahort-aeal"  solutions  for  incompressible  aeala,  which 
account  for  circumferential  fluid  flow  due  to  wall  ahear  atreaaes  but 
not  preaaure  perturbationa.  The  analyaia  employa  a  bulk-flow  assumption 
and  aooounta  for  fluid  prerotation  aa  ia  enters  the  aeal.  Childs'  [7] 
incompressible  aeal  analyaia  provides  "finite-length  "  solutions,  in 
whioh  both  ahear  and  preaaure- Induced  flow  are  included.  Childs' 
utilizes  Hira'  [8]  turbulent  bulk-flow  model,  and  aooounta  for  inlet 
swirl  aa  well  a a  perturbationa  in  axial  and  oiroumferentlal  Reynolds 
numbers  due  to  clearance  perturbationa. 

Compressible  flow  in  seals  haa  been  analyzed  by  Fleming  [9,  10]  and 
Nelson[11,  12].  Fleming  presents  a  short  aeal  solution  for  the  leakage, 
direot  stiffness,  and  direct  damping  coefficients  for  straight  and 
tapered,  smooth,  annular  gas  seals,  but  does  not  include  the 
oroaa-coupled  damping  terms.  Nelson,  whose  analysis  is  used  for 
comparison  in  this  report,  analyzes  both  smooth  and  surfaoe-roughened 
annular  seals  in  the  straight  and  tapered  configurations.  An  outline  of 
Nelson's  analyaia  is  included  in  the  seotion  that  follows. 


NELSON'S  ANALYSIS 


Nelson  [11,  12]  has  developed  an  analysis  whioh  provides  both 

static  and  dynamic  results  for  annular  gas  seals.  The  statio  results 
Include  fluid  leakage  through  the  seal,  pressure  gradient  along  the  seal 
axis,  and  the  fluid  axial  and  circumferential  velocities  through  the 
seal.  Dynamic  data  provided  by  the  analysis  consists  of  the 
rotordynamio  ooeffioients  (direot  and  cross-coupled  stiffness  and 
damping  terms)  for  small  rotor  motion  about  a  centered  position 
(equation(2) ) .  Nelson  assumes  that  the  added-inass  terms  are  negligible 
for  gas  seals,  and,  hence,  equation^)  is  written 

,/ 


K 

-k 


c 

“C 


(3) 


Nelson  utilizes  a  modified  Kirs'  [8]  turbulent  bulk-flow  fluid 
model  to  develop  governing  axial  and  circumferential  momentum  equations, 
and  his  model  is  completed  by  the  continuity  and  energy  equations. 
Kirs'  model  defines  the  wall  shear  stress  tw  as 

xw  -  1/2  pUm2  noUpUnjH  /  u)"10  -  1/2  pUm2  noRamo  (4) 
where  Um  is  the  mean  flow  velocity  relative  to  the  surface  upon  which 
the  shear  stress  acts,  and  H  is  the  local  seal  clearance.  Kirs' 
formulation  assumes  that  the  surface  roughness  is  the  same  on  the  rotor 
and  stator.  However,  if  the  bulk-flow  velocities  relative  tu  the  rotor 
and  stator  are  substituted  in  equation  (4),  the  shear  stresses  at  the 
rotor  and  stator  are,  respectively, 
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(5) 


.  tr  -  1/2  pUr2  nr(2pUrfl  /  u)mr 
ts  -  1/2  pUs2  ns(2pUaH  /  u)ms 
Hence,  different  surface  roughnesses  in  the  seal  elements  can  be 
accounted  for  via  the  empirical  coefficients  mr,  nr  and  ms, ns  for  the 
rotor  and  stator  surfaces.  These  coefficients  may  be  calculated  from 
statio-pressure-gradient  teat  data,  and  are  then  provided  as  input 
parameters  for  Nelson's  analysis. 

Assuming  small  motion  of  the  rotor  about  a  centered  position, 
Nelson  uses  a  perturbation  analysis  similar  to  that  employed  by  Childs 
[7]  to  develop  zeroth  and  first-order  perturbation  equations.  The 
zoroth-order  solution  represents  a  zero-eccentricity  flow  condition, 
with  rotor  rotation  but  without  precession.  This  solution  is  iterative 
and  yields  the  mass-leakage  flow  rate,  and  the  axial  distribution  of 
pressure,  axial  velocity,  density,  and  circumferential  velocity. 

An  iterative  solution  saheme  is  employed,  using  initial  guesses  for 
the  zeroth-order  seal  entrance  Mach  number  and  entrance  pressure-loss 
coefficient.  The  entrance-loss  relationship  is  defined  by 

Y/(Y-1) 

po(0)  -  1  /  (1  +  C  (Y-l ) (k+1 )Mo2(0) ]  /  2}  (6) 


where  po(0)  is  the  seal  entrance/reservoir  pressure  ratio  and  Mo(0)  is 
the  entrance  Maoh  number.  The  entrance  Maoh  number  is  iteratively 
adjusted,  and  the  loss  coefficient  1<  is  recalculated  according  to  a 
curve  fit  by  Deissler  [13] 


E  +  1  -V5. 3  /  logio  «a  (?) 
which  is  plotted  in  Fig.  3.  At  axial  Reynolds  numbers  above  200,000,  k 
is  equated  to  zero.  The  iterative  solution  procedure  for  Mo(0)  and  k 
continues  until  either: 


1)  the  Mach  number  at  the  exit  reaches  unity  and  the  exit  pressure 
is  greater  than  the  sump  pressure  (choked  flow),  or 

2)  the  exit  pressure  equals  the  sump  pressure  and  the  exit  Mach 
number  Is  less  than  unity  (unohoked  flow). 

The  pressure,  density,  and  velocity  distribution  and  their 
derivatives  whioh  are  determined  in  the  zeroth-order  solution  are  used 
in  defining  coefficients  of  the  first-order  perturbation  equations. 
These  equations  define  the  pressure,  density,  and  axial  and 
circumferential  velooity  perturbations  due  to  rotor  motion,  and  are 
transformed  to  sixteen  ordinary  differential  equations  in  the  axial 
coordinate  z.  The  four  physical  boundary  conditions  required  for 
the  solution  of  these  aquations  depend  on  the  perturbation  conditions 
that  are  specified  at  the  3eal  entrance  and  exit.  The  inlet 
circumferential  velooity  perturbation  is  zero.  Expansion  or  the 
entrance  pressure-loss  relationship  of  equation  (6)  yields  a  second 
boundary  condition.  For  ohoked  flow,  the  first-order  perturbation  in 
the  exit  Maoh  number  is  zero,  while  for  unchoked  flow,  the  firat-order 
perturbation  in  the  exit  pressure  is  zero. 

Application  of  these  boundary  conditions  and  numerical  integration 
of  the  ordinary  differential  equations  provides  the  first-order 
solution.  Integration  of  the  first-order  pressure  solution  along  and 
around  the  seal  periphery  yields  the  direct  and  croes-coupled  stiffness 
and  damping  coefficients,  K,  k,  C,  and  c,  respectively. 

The  input  parameters  whioh  can  be  varied  in  Nelson's  analysis 
incl udes 

1)  reservoir  pressure  and  temperature, 

2)  sump  pressure, 


3)  seal  geometry  (i.e.  radius,  length,  clearances), 

4)  rotor  rotational  speed  and  precession  rate, 

5)  entrance  circumferential  velocity  of  fluid, 

6)  rotor  and  stator  surface  roughness  (Hirs  constants), 

7)  empirical  entranoe-loss  relationship,  e.g.,  Deissler's, 
equation  (16),  and 

8)  fluid  viscosity,  gas  constant,  and  ratio  of  speoifio  heats. 

Zt  is  apparent  that  a  large  amount  of  theoretical  data  can  be 
generated  to  determine  the  influence  that  these  various  parameters  have 
on  the  fluid  foroes  in  annular  gas  seals.  However,  there  is  a  lack  of 
experimental  data  with  which  to  oompare  the  results  of  Nelson's 
analysis.  Currently,  test  results  due  to  Wachter  and  Benokert  [14]  exist 
for' labyrinth  seals,  a  speoial  class  of  non-eontaoting  seals  which  have 
stepped  surfaces  or  "teeth"  on  the  rotor,  stator,  or  both.  Experimental 
results  for  smooth  and/or  surface-roughened  gas  seals  are  limited  to 
data  for  honeyoomb  seals  also  published  by  Wachter  and  Benokert.  Hence, 
the  need  for  a  teat  apparatus  which  can  be  used  to  study  the  effects  of 
the  same  variables  provided  for  in  Nelson's  analysis  is  obvious.  The 
experimental  data  generated  by  such  an  apparatus  would  be  valuable  for 
comparison  to  both  Nelson’s  theories  and  others  which  may  be  developed 


in  the  future. 
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TEST  CONCEPTS 

A  number  of  test  programs  have  been  implemented  to  measure  the- 
stabilizing  and  destabilizing  fluid  forces  which  are  developed  by 
turboraachinery  elements.  Some  are  concerned  mainly  with  the  study  of 
seal  forces,  while  others  examine  the  forces  developed  by  centrifugal 
pump  impellers.  In  eaoh  case,  reaction  force  and  relative  motion 
measurements  are  used  for  rotordynamio  coefficient  identification.  Four 
general  approaches  have  been  employed,  and  will  be  reviewed  here. 

Waohter  and  Benckert  [14]  employ  a  statio  displacement  method  for 
determining  stiffness  coefficients.  In  this  method,  as  shown  in  Fig.  4, 
the  rotor  is  displaced  statioally  to  some  measured  ecoentrio  position 
while  a  pressure  differential  forces  the  working  riuid  past  the  seal. 
By  measuring  the  reaction  force  components  whioh  are  parallel  and 
perpendioular  to  the  static  displacement  vector,  the  direct  and 
cross-coupled  stiffnesses  can  be  determined.  Referring  to  equation  (2) 
for  small  rotor  motion  about  a  centered  position,  a  statio  rotor 
displacement  in  the  X-direction  yields 

K  -  -Fx/e0  -  k  -  (8) 
Since  this  static  displacement  method  has  no  dynamic  motion,  no  damping 
or  added-mass  terms  can  be  evaluated. 

A  second  approach  to  rotordynamio  coefficient  identification  is 
utilized  by  Childs[15].  Depicted  in  Fig. 5,  this  method  uses  a  circular 
orbit  of  the  rotor  within  the  seal.  The  rotor  is  mounted  eccentrically 
on  a  shaft  whioh  rotates.  Thus,  the  rotor  precesses  in  a  circular  orbit 
at  the  same  rate  and  direction  as  shaft  rotation.  •  This  synchronous 
precession  provides  for  the  determination  of  the  radial  and  tangential 
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components  of  the  seal  reaction  force.  The  Fr  and  F  components  are 
obtained  through  integration  of  the  measured  pressure  distribution  along 
and  around  the  seal  periphery.  Expressing  measured  rotor  motion  as 

X  -  e0  oos(uit) 

(9) 

Y  «  e0  sin(uit) 

for  small  oiroular  orbit  of  radius  e0  and  preoessional  frequency  u-ft, 
and  substituting  into  equation  (2)  yields  the  radial  and  tangential 
foroe  ooefflooient  definitions 

* 

Fr  /  e0  -  Mu1  -  ooj  -  K  ■  -Kef  +  Mgf.u)* 

(10) 

F  /  Sq  *  k  —  Cu  “  "Cgf .w 

where  the  oross-ooupled  mass  ooeffioient  is  assumed  negligible  with 
respeot  to  the  influence  of  k  and  C.  Beoause  the  oross-ooupled 
coefficients  k  and  o  are  linear  funotions  of  oj,  identification  of  the 
individual  dynamio  ooeffioients  is  not  possible  in  this  method.  However, 
equivalent  direct  stiffness,  damping,  and  added-mass  ooeffioients  oan  be 
oaloulated  as  indioated  in  equation  (10.). 

Independent  rotation  and  precession  control ,  as  shown  in  Fig.  6,  is 
a  third  testing  method  which  is  currently  employed  both  in  impeller  and 
seal  studies  [16],  [17],  [18].  Various  means  are  used  to  produce  a 

oiroular  orbit  (precession)  of  the  rotor  or  impeller  at  a  rate  different 
from  its  rotational  speed.  For  a  small  circular  orbit  of  radius  e0  and 
preoessional  frequency  fl  ,  the  measured  preoessional  motion  of  the  rotor 
is 

X  -  e0  oos(nt) 

Y  -  e0  sin(Pt) 

The  Fy  and  Fy  reaction  force  components  are  measured  and  oan  be 


components  of  the  seal  reaction  force.  The  Fr  and  F  components  are 
obtained  through  Integration  of  the  measured  pressure  distribution  along 
and  around  the  seal  periphery.  Expressing  measured  rotor  motion  as 

X  -  e0  oos(wt) 

(9) 

X  -  e0  sin(wt) 

for  small  oiroular  orbit  of  radius  e0  and  preoessional  frequency 
and  substituting  into  equation  (2)  yields  the  radial  and  tangential 
force  coefficoient  definitions 

* 

Ff.  /  e0  ■  Mu2  -  cw  -  K  *  -Kef  +  Mef.(D2 


(10) 


F  /  e0  ■  k  -  Cu  ■  -C8f.w 

where  the  cross-coupled  mass  coefficient  is  assumed  negligible  with 
respeot  to  the  influence  of  k  and  C.  Because  the  oross-ooupled 
coefficients  k  and  o  are  linear  functions  of  u,  identification  of  the 
individual  dynamia  coefficients  Is  not  possible  in  this  method.  However, 
equivalent  direct  stiffness,  damping,  and  added-mass  coefficients  can  be 
calculated  as  indicated  in  equation  (10.). 

Independent  rotation  and  precession  control,  as  3hown  in  Fig.  6,  is 
a  third  testing  method  which  is  currently  employed  both  in  impeller  and 
seal  studies  C 1 6] ,  [17],  [18].  Various  means  are  used  to  produoe  a 
oiroular  orbit  (preoession)  of  the  rotor  or  impeller  at  a  rate  different 
from  its  rotational  speed.  For  a  small  circular  orbit  of  radius  e0  and 
preoessional  frequency  fl  ,  the  measured  preoessional  motion  of  the  rotor 
is 


X  -  e0  cos(nt) 

X  -  e0  sin(at) 

The  Fy  and  Fy  reaction  foroe  components  are  measured  and  can  be 


expressed 


FX  *  sin(flt)  +  FXq  cos  (Jit) 
Fy  •  Fys  sin(flt)  +  Fyc  oos(JJt) 


By  substituting  these  expressions  into  equation  (2)  and  equating 


ooeffioients  of  sine  and  oosine  terms,  the  following  equations  are 


obtained 


-  FXc  /  e0  ■  K  +  oJl  -  Mil* 


Fyc  /  «o  “  +  cn  *mnS 


-  Fxs  /  e0  -  k  -  cn  -  mJla 


FyS  /  e0  -  K  +  cQ  +  MS2J 


Henoe,  by  measuring  the  reaotion  force  components  and  rotor  motion  at 


two  different  precession  frequencies,  eight  equations  in  six  unknowns 


are  obtained,  and  the  rotordynamio  ooeffioients  oan  be  calculated. 


A  fourth  testing  method  has  been  used  by  lino  and  Kaneko  [19]  for 


determining  dynamio  ooeffioients,  and  this  same  method  is  employed  at 


the  TAMU  gas  seal  test  facility.  An  external  hydraulio  shaker  is  used 


to  Impart  translatory  harmonic  motion  to  the  rotating  seal,  and  rotor 


motion  relative  to  the  stator  and  the  reaction  force  components  acting 


on  the  stator  are  measured. 


Fig.  7  shows  the  manner  in  which  the  rotor  could  be  positioned  and 


osoillated  in  order  to  identify  the  dynamic  ooeffioients  of  the  seal  for 


small  motion  about  e0.  If  the  added-mass  terms  are  assumed  negligible, 


equation  (1)  is  rewritten 


Fxl  KXX<eo)  Kxy(e0)|  \  X  /  Cxx(e0)  Cxy(eo)  I  1 X 


KYX<eo>  KYY <e0)  |  |y|  |CYx(eo)  CYYUo)  I  I* 


XX'eO'  XY'eo 


'YX'eo'  YY ' eo 


First,  harmonic  horizontal  motion  of  the  rotor  is  assumed,  where 
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X  -  e0  +  A  sin(flt)  +  B  ooa(at) 

X  -  Afl  ooa(flt)  -  BQ  ain(flt) 

Y  •  Y  ■  0 

Thia  yields  amall  motion  parallel  to  the  atatic  eccentricity  vector, 
where  a  ia  the  shaking  frequency.  In  a  aimilar  faahion,  the  X  and 
Y-direotion  force  oomponenta  oan  be  expreaaed 

Fv  -  Fvc  sin(flt)  +  Fyq  ooa(at) 

(HO 

Fy  -  Fys  ain(at)  +  Fyc  ooa(at) 

Subatituting  theae  expreasiona  into  equation  (13)  and  equating 
ooeffioienta  of  aine  and  cosine  terma  yielda  the  following  four 
equations 

FXS“  KXX  A  “  CXX  B 
FXC  *  KXX  B  +  CXX  A 

(15) 

FyS  -  KyX  A  -  CyX  B 
FyC  -  Kyx  B  +•  Cyy  A 

Solving  thia  ayatera  of  four  equations  in  four  unknowns  defines  the 
dynamic  ooeffioienta  aa 

KXx<e0)  -  (FXC  B  +  Fxs  A)  /  (A*  +  B2) 

Kyx(Eo)  ■  (Fys  A  +  Fyc  B)  /  (A1  +  B2 ) 

(16) 

CXx<e0)  -  (FXC  A  -  Fxs  B)  /  n(Aa  ♦  Bl) 

Cyx(Co)  "  <FYC  A  “  FyS  B)  /  fl(  A2  ♦  B2) 

Therefore,  by  measuring  the  reaction  forces  due  to  known  rotor 
motion,  determining  the  Fourier  coefficients  (A,  B,  FXs,  FXc»  Fys,  Fyc)> 
and  substituting  into  the  above  definitions,  the  indicated  dynamic 
ooeffioienta  can  be  identified.  If  the  rotor  ia  shaken  about  a  centered 
position,  then  the  process  is  complete.  Since  the  linearized  model  has 
akew-symmetric  stiffness  and  damping  matrices,  all  of  the  coefficients 


are  idantifled.  It,  however,  the  rotor  is  shaken  about  an  eooentrio 
position  as  initially  postulated,  then  it  must  bo  shaken  vertically 
about  that  sane  point  in  order  to  oomplete  the  identification  process. 
Assuming  harmonic  vertloal  motion  of  the  rotor,  as  defined  by 

X  -  e0,  X  -  0, 

Y  ■  A  sin(flt)  +  B  cos(Ot),  and 

Y  -  AQ  cos(flt)  -  BO  sinfnt), 

yields  osoillatory  motion  that  is  perpendicular  to  the  assumed  static 
eooentricity  veotor.  A  similar  process  as  before  results  in  the 
coefficient  definitions 

KYY<«o)  -  <FXS  A  +  FXC  B)  /  (A*  +  B4) 

Kxy<c0)  -  "(Fyc  B  ♦  Fys  A)  /  (A*  ♦  B4) 

(17) 

C*y<«o)  "  <FXC  A  -  Fxs  B)  /  Q(  A4  ♦  B4) 
cXY^eo)  “  (pYS  0  "  Fyc  A)  f  B(A4  + 


All  eight  dynamio  coefficients  are  thus  determined  by  alternately 
shaking  the  rotor  at  one  frequency  fl  in  directions  which  are  parallel 
and  perpendicular  to  the  static  eccentricity  vector. 


TEST  APPARATUS  OVERVIEW 


Detailed  design  of  the  TAMU  gas  seal  apparatus  was  carried  out  by 
J.B.  Dressman  of  the  University  of  Louisville.  It  is  of  the  external 
shaker  configuration,  and  the  dynamic  coefficient  identification  process 
is  as  described  in  the  latter  part  of  the  preceding  seotion. 

Considering  both  the  coefficient  identification  process  and 
Nelson's  analysis,  some  objectives  for  the  design  of  the  test  apparatus 
are  apparent.  First,  in  order  to  determine  the  dynamic  coefficients, 
the  apparatus  must  provide  for  the  neoessary  rotor  motion  within  the 
seal,  and  measurement  of  the  reaot ion- force  components  due  to  this 
motion  must  be  possible.  Secondly,  it  would  be  advantageous  (for 
purposes  cf  comparison)  if  the  apparatus  oould  provide  the  same  variable 
seal  parameters  afforded,  by  Nelson's  analysis  (i.e.,  pressures,  seal 
geometry,  rotor  rotational  speed,  fluid  prerotation,  and  rotor/stator 
surface  roughness).  With  this  capability,  the  influenoe  of  each 
independent  parameter  could  be  examined  and  compared  for  correlation 
between  theoretical  predictions  and  experimental  results. 

With  these  design  objectives  in  mind,  the  disoussion  of  the  test 
apparatus  is  presented  in  three  sections.  The  first  section,  Test 
Hardware,  describes  how  the  various  seal  parameters  are  physically 
executed  and  controlled.  For  example,  the  manner  in  which  the  dynamic 
"shaking"  motion  of  the  seal  rotor  is  achieved  and  controlled  is 
described  in  this  section.  The  second  section,  Instrumentation, 
describes  how  these  controlled  parameters,  such  as  rotor  motion,  are 
measured.  Finally,  the  Data  Acquisition  and  Reduction  section  explains 
how  these  measurements  are  used  to  provide  the  desired  information. 
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TEST  HARDWARE 


This  section  deals  only  with  the  mechanical  components  and 
operation  of  the  test  apparatus.  It  is  intended  to  provide  answers  to 
the  following  questions: 

1)  How  is  the  static  position  of  the  seal  rotor  controlled? 

2)  How  is  the  dynamic  motion  of  the  rotor  executed  and  controlled? 

3)  How  is  compressed  air  obtained  and  supplied  to  the  apparatus, 
and  how  is  the  pressure  ratio  across  the  seal  controlled? 

4)  How  is  the  incoming  air  prerotated  before  it  enters  the  seal? 

5)  How  are  the  seal  rotor  and  stator  mounted  and  replaoed? 

6)  How  is  the  seal  rotor  driven  (rotated)? 

Recalling  the  rotordynamio  ooeffioient  identification  process 
described  earlier,  the  external  shaker  method  requires  that  the  seal 


rotor  be  set  in  some  static  position  and  then  be  oscillated  about  that 
point.  The  teat  apparatus  meets  those  requirements  by  providing 
Independent  statio  and  dynamic  displacement  oontrol,  which  are  desoribed 
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below. 

Static  Displacement  Control.  The  test  apparatus  is  designed  to  provide 
oontrol  over  the  static  eccentricity  position  both  horizontally  and 
vertically  within  the  seal.  The  rotor  shaft  is  suspended 
pendulum-fashion  from  an  upper,  rigidly  mounted  pivot  shaft,  as  shown  in 
Figs.  8  and  9.  This  arrangement  allows  a  side-to-side  (horizontal) 
motion  of  the  rotor,  and  a  oam  within  the  pivot  shaft  allows  vertical 
positioning  of  the  rotor. 

The  cam  which  controls  the  vertioal  position  of  the  rotor  Is  driven 
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by  a  remotely-operated  DC  gearhead  motor,  allowing  aocurate  positioning 
of  the  rotor  during  testing.  Horizontal  positioning  of  the  rotor  is 
accomplished  by  a  Zonio  hydraulic  shaker  head  and  master  controller, 
which  provide  Independent  statio  and  dynaoio  displacement  or  force 
oontrol.  The  shaker  head  is  mounted  on  an  I-beam  support  structure,  and 
oan  supply  up  to  4450  N  (1000  lbf)  static  and  4450  N  dynamic  force.  As 
illustrated  in  Fig.  8,  the  shaker  head  output  shaft  aots  on  the  rotor 
shaft  bearing  housing,  and  works  against  a  return  spring  mounted  on  the 
opposite  side  of  the  bearing  housing.  The  return  spring  maintains 
oontact  between  the  shaker  head  shaft  and  the  bearing  housing,  thereby 
preventing  hammering  of  the  shaker  shaft  and  the  resulting  loss  of 
oontrol  over  the  horizontal  motion  of  the  rotor. 

Dynamio  Diaplaoement  Control.  The  dynamic  motion  of  the  seal  rotor 
within  the  stator  is  horizontal.  In  addition  to  controlling  the  statio 
horizontal  position  of  the  rotor,  the  Zonic  shaker  head  moves  the  rotor 
through  horizontal  harmonia  oscillations  as  the  test  is  run.  A  Wavetek 
funotion  generator  provides  the  sinusoidal  input  signal  to  the  Zonic 
controller,  and  both  the  amplitude  and  frequency  of  the  rotor 
oscillations  are  controlled. 

Although  the  test  rig  design  provides  for  dynamic  motion  of  the 
rotor  only  in  the  horizontal  X-direotion,  all  of  the  coefficients  for 
either  seal  model  (equation  (3)  or  (13))  can  still  be  determined.  As 
Fig.  10  shows,  the  required  rotor  motion  perpendicular  to  the  static 
ecoentricity  vector  can  be  accomplished  in  an  equivalent  manner  by 
statically  displacing  it  the  same  amount  (e0)  in  the  vertical  direction 
and  continuing  to  shake  horizontally. 

In  addition  to  providing  control  over  the  rotor's  static  position 


and  dynamio  motion,  the  test  apparatus  allows  other  seal  parameters  to 
be  controlled  independently,  providing  insight  into  the  influence  these 
parameters  have  on  seal  behavior.  These  parameters  coincide*  with  the 
variable  input  parameters  for  Nelson's  analysis,  and  they  include: 

1 )  pressure  ratio  across  the  seal , 

2)  prerotsMoh  of  the  incoming  fluid, 

3)  seal  corn  i duration,  and 

4)  rotor  rotational  speed. 

Pressure  Ratio.  The  inlet  air  pressure  and  attendant  mass  flow  rate 
through  the  seal  are  controlled  by  an  eleatrlcover-pneumatloally 
actuated  Masoneilan  Camflex  II  flow  oontrol  valve  located  upstream  of 
the  test  seotlon.  An  Ingersoll-Rand  SSR-2000  single  stage  sorew 
compressor  rated  at  34  m'/min  9  929  kPa  (1200  sofm  6  120  psig)  provides 
oorapresaed  air,  which  is  then  filtered  and  dried  before  entering  a  surge 
tank.  Losses  through  the  dryers,  filters,  and  piping  result  In  an 
aatual  maximum  inlet  pressure  to  the  test  section  of  approximately 
722  kPa  (90  psig)  and  a  maximum  flow  rate  of  27  m’/min  (950  sofm).  A 
four-inoh  inlet  pipo  from  the  surge  tank  supplies  the  test  rig,  and 
after  passing  through  the  seal,  the  air  exhausts  to  atmosphere  through  a 
manifold  with  muffler. 

Inlet  Clroumf erentlal  Velocity  Control.  In  order  to  determine  the 
effects  of  fluid  rotation  on  the  rotordynamio  coefficients,  the  test  rig 
design  also  allows  for  prerotation  of  the  Incoming  air  as  it  enters  the 
seal.  This  prerotation  introduces  a  oiroumferentlal  component  to  the 
air  flow  direotion,  and  is  accomplished  by  guide  vanes  which  direct  and 
acoelerate  the  flow  towards  the  annulus  of  the  seal.  The  vanes 
are  maohined  from  brass  disku,  and  Fig.  11  illustrates  the  vane 


30 


oonf iguration.  Three  seta  of  guide  vanes  are  available!  one  rotates  the 
flow  in  the  direotion  of  rotor  rotation,  another  introduces  no  flu^d 
rotation,  and  the  third  rotates  the  flow  opposite  the  direotion  of  rotor 
rotation. 

Seal  Configuration.  The  design  of  the  test  rig  permits  the  installation 
of  various  rotor/stator  combinations.  As  shown  in  Figs. 12-15,  the 
stator  is  supported  in  the  test  seotion  housing  by  three  Kistler  quartz 
load  cells  in  a  trihedral  oonf iguration.  Figs. 12  and  13  show  the 
amooth-rotor/smooth-stator  seal,  while  the  smooth-rotor/honeyoomb-stator 
seal  is  illustrated  in  Figs. 14  and  15.  The  seal  rotor  is  press-fitted 
and  seoured  axially  by  a  bolt  oirole  to  the  rotor  shaft.  Seals  with 
different  geometries  (i.e.,  clearanoea,  tapers,  lengths)  oan  be  tested, 
as  well  as  seals  with  different  surfaoe  roughnesses.  The  replacement  of 
these  rotor/stator  combinations  oan  be  accomplished  with  minimal 
downtime. 

Rotational  Speed.  A  Westinghouse  50-hp  variable-speed  eleotrio  motor 
drives  the  rotor  shaft  through  a  belt-driven  Jaokshaft  arrangement. 
This  shaft  is  supported  by  two  sets  of  Torrington  hollow-roller  bearings 
[20].  These  bearings  are  extremely  preoise,  radially  preloaded,  and 
have  a  predictable  and  repeatable  radial  stiffness.  Axial  thrust  due  to 
the  pressure  differential  across  the  seal  is  absorbed  by  a  flat, 
roller-type,  caged  thrust  bearing  at  the  rear  of  the  rotor.  Both  the 
shaft  and  thrust  bearings  are  lubricated  by  a  positive-displaoement 
gear-type  oil  pump. 

Different  Jaokshaft  drive-pulleys  oan  be  fitted  to  provide  up  to  a 
4i1  speed  increase  from  motor  to  rotor  shaft,  which  would  result  Ln 
a  rotor  shaft  speed  range  of  0-21,200  cpm.  Current  design  limitations, 
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however,  prevent  the  attainment  of  this  upper  rotational  speed.  High 
bearing  temperatures,  reduotion  of  interference  in  the  rotor-rotor  shaft 
fitment  due  to  inertia-induoed  radial  growth  of  the  rotor  inside 
diameter,  and  exoessive  stresses  in  the  drive-pulleys  have  served  to 
limit  shaft  speed.  The  highest  rotational  speed  attained  at  the  time  of 
this  writing  is  8500  opm,  although  design  modifications  to  allow  higher 
speeds  are  under  investigation. 

To  oonolude  this  disoussion  of  the  test  hardware,  two  views  of  the 
complete  test  apparatus  are  inoluded.  Fig.  16  shows  the  assembled  rig, 
while  an  exploded  view  is  provided  in  Fig.  17. 
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INSTRUMENTATION 

Having  discussed  what  aeal  parameters  can  be  varied,  and  how  the 
variations  are  implemented,  the  measurement  of  their  respective  effeots 
oan  now  be  daaoribed.  The  types  of  measurements  which  are  made  oan  be 
grouped  into  three  categories: 

1 )  rotor  motion, 

2)  reaotion-foroe  measurements,  and 

3)  fluid  flow  measurements. 

These  categories  are  desoribed  individually  in  the  sections  that  follow. 
Rotor  Motion  Measurements.  The  position  of  the  sea!)  rotor  within  the 
stator  is  monitored  by  two  Bently-Nevada  eddy-current  proximity  probes, 
mounted  in  the  teat  seotion  housing.  These  probes  are  located  90  degrees 
apart,  and  correspond  to  the  X  and  Y-  directions.  The  proximity  probes 
are  used  to  determine  the  statlo  position  and  dynamio  motion  of  the 
rotor,  and  their  resolution  is  0.0025  mm  (0.1  mil). 

Reaction-Force  Measurements.  Reaction  forces  arise  due  to  the  static 
position  and  dynamio  motion  of  the  seal  rotor  within  the  stator.  The 
reaotion  forota  (Fy,  Fy)  exerted  on  the  stator  are  measured  by  the  three 
Kistler  quartz  load  cells  whioh  support  the  stator  in  the  test  seotion 
housing.  When  the  rotor  is  shaken,  vibration  is  transmitted  to  the  test 
seotion  housing,  both  through  the  thrust  bearing  and  through  the  housing 
mounts.  The  acceleration  of  the  housing  and  stator  generates  unwanted 
Inertial  "ma"  foroes  whioh  are  sensed  by  the  load  cells,  in  addition  to 
those  pressure  foroes  developed  by  the  relative  motion  of  the  seal  rotor 
and  stator.  For  this  reason,  PCB  piezoelectric  accelerometers  with 
integral  amplifiers  are  mounted  in  the  X  and  Y-direotions  on  the  stator, 
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as  shown  in  Figs.  13  and  15.  These  aooels  allow  a  (stator  mass)  x 
(stator  aooeleration)  subtraction  to  the  forces  (Fx*  Fy)  indicated  by 
the  load  cells.  With  this  oorreotion,  which  is  described  more  fully  in 
the  next  section,  only  the  pressure  foroes  due  to  relative  seal  motion 

are  measured. 

Force  measurement  resolution  is  a  function  of  the  stator  mass  and 
the  resolution  of  the  load  cells  and  accelerometers.  Aocelerometer 
resolution  is  0.005  g,  which  must  be  multiplied  by  the  stator  mass  in 
order  to  obtain  an  equivalent  foroe  resolution.  The  masses  of  the 
stators  used  in  the  test  program  reported  hero  are  11.4  kg (25. 2  lb)  and 
3.94  kg (8. 69  lb),  corresponding  to  the  smooth  and  honeycomb  stators, 
respectively.  Henoe,  foroe  resolution  for  the  aooelerometers  is  0.560  N 
(0.126  lb)  and  0.191  N  (0.043  lb),  for  each  stator,  respectively. 
Resolution  of  the  load  cells  is  0.089  N(0.02  lb).  Therefore,  the 
resolution  of  the  foroe  measurement  is  limited  by  the  aooelerometers . 
With  a  stator  with  leas  mass,  and/or  aooelerometers  with  greater 
sensitivity,  foroe  resolution  oouid  be  improved. 

Fluid  Flow  Measurements.  Fluid  flow  measurements  include  the  leakage 
(mass  flow  rate)  of  air  through  the  seal,  the  pressure  gradient  along 
the  seal  axis,  the  inlet  fluid  circumferential  velocity,  and  the 
entranoe  pressure  loss. 

Leakage  is  measured  with  a  Fischer  &  Porter  vortex  flowmeter 
looated  in  the  piping  upstream  of  the  test  section.  Resolution  of  the 
flowmeter  is  0.0014  m1  (0.05  acf),  and  pressures  and  temperatures  up  and 
downstream  of  the  meter  are  measurod  for  mass  flow  rate  determination. 

For  measurement  of  the  axial  pressure  gradient,  the  stator  has 
pressure  taps  drilled  along  the  length  of  the  seal  in  the  axial 


direction.  These  pressures,  as  well  as  all  others,  are  measured  with 
a  0-1.034  MPa  (0-150  psig)  Soani valve  differential-type  pressure 
transducer  through  a  48  port,  remotely-controlled  Soani valve  model  J 
soanner.  Transducer  resolution  is  0.552  kPa  (0.08  psi). 

In  order  to  determine  the  oiroumferential  velocity  of  the  air  as 
it  enters  the  seal,  the  static  pressure  at  the  guide  vane  exit  is 
measured.  This  pressure,  in  conjunction  with  the  measured  flowrate  and 
inlet  air  temperature,  is  used  to  calculate  a  guide  vane  exit  Mach 
number.  A  oompressible  flow  continuity  equation 

m  -  Pax  Aex  M0x  C(T/RgTt)  (1  *  (Y-1)Max2  /  2)]  1 '2  (18) 
is  rearranged  to  provide  a  quadratic  equation  for  Mex 

Mex2  -  {-1  +  1  ♦  4((Y-1)/2Y)  (m  RgTt  /  pex  Asx)2  }  /  (7-1)  (19) 
where  Y  1s  the  ratio  of  speoific  heats  and  Rg  is  the  gas  constant  for 
air,  Tt  is  the  stagnation  temperature  of  the  air,  pex  is  the  static 
pressure  at  the  vane  exit,  and  Aex  is  the  total  exit  area  of  the  guide 
vanes.  Since  all  of  the  variables  in  the  equation  are  either  known  or 
measured,  the  vane  exit  Mach  number,  and  therefore  the  velooity,  can  be 
found. 

In  order  to  determine  the  oiroumferential  component  of  this  inlet 
velooity,  a  flow  turning  angle  correction,  in  accordance  with  Cohen[21], 
is  employed.  The  correction  has  been  developed  from  guide  vane  cascade 
tests,  and  accounts  for  the  fact  that  the  fluid  generally  is  not  turned 
through  the  full  angle  provided  by  the  shape  of  the  guide  vanes.  With 
this  flow  deviation  angle  calculation,  the  actual  flow  direction  of  the 
air  leaving  the  vanes  (and  entering  the  seal)  can  be  determined.  Hence, 
the  magnitude  and  direction  of  the  inlet  velocity  is  known,  and  the 
appropriate  component  is  the  measured  inlet  circumferential  velocity. 


The  entrance  preaaure-loss  coefficient,  defined  in  equation  (6),  is 
determined  from  the  measured  pressures  just  upstream  of  and  Just  inside 
the  seal.  An  entranoe  Mach  number  is  calculated  in  the  same  manner  as 
outlined  previously,  using  the  measured  pressure  immediately  inside  the 
seal  and  the  annular  area  between  the  rotor  and  stator.  This  entranoe 
Maoh  number,  and  the  ratio  of  the  seal  entrance/ guide  vane  exit 
pressures  are  substituted  into  equation  (6),  and  the  entrance  loss 
coefficient,  k,  is  determined. 


DATA  ACQUISITION  AND  REDUCTION 


With  the  preceding  explanations  of  how  the  seal  parameters  are 
varied,  and  how  these  parameters  are  measured,  the  discussion  of  how  the 
raw  data  is  processed  and  Implemented  oan  begin.  Data  acquisition  is 
dlreoted  from  a  Hewlett-Packard  901 6  (16-bit)  computer  with  disk  drive 
and  9.8  megabyte  hard  disk.  The  computer  controls  an  H-P  6940B 
multi  programmer  whioh  has  12-bit  A/D  and  D/A  converter  boards  and 
transfers  oontrol  commands  to  and  teat  data  from  the  instrumentation. 

As  was  previously  stated,  the  major  data  groups  are  seal 
motlon/reaotion  force  data  and  fluid  flow  data.  The  motion/reaotion 
force  data  are  used  for  dynamlo  coefficient  identification.  The 
hardware  involved  lnoludes  the  load  cells,  aooelerometers,  X-direotlon 
motion  probe,  a  Sensoteo  analog  filter  unit,  a  tunable  bandpass  niter, 
and  the  A/D  converter.  The  operation  of  these  components  is  illustrated 
in  Fig. 18,  and  their  outputs  are  used  in  a  serial  sampling  scheme  whioh 
provides  the  oomputer  with  the  desired  data  for  reduction.  Rooalling 
the  discussion  of  the  reaction  force  measurements  in  the  preooding 

i 

seotion,  a  (stator  mass)  x  (stator  acceleration)  subtraction 
from  the  indicated  load  aell  forces  is  necessitated  due  to  vibration  of 
the  stator  and  test  seotion  housing.  This  subtraction  is  performed 
with  an  analog  clroult,  and  results  in  corrected  Fx  and  Fy  force 
components  due  to  relative  seal  motion. 

The  foroed  osoillotory  shaking  motion  of  the  seal  rotor  is  the  key 
to  the  operation  of  the  serial  synchronous  sampling  (SSS)  routine  which 
is  employed.  The  frequency  of  the  rotor  osoillation  is  set  by  a  funotion 
generator,  and  rotor  motion  is  sensed  by  the  X, -direction  motion  probe. 
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Th«  notion  signal  is  flltsrcd  by  ths  narrow  bandpass  filter,  and  is  used 
as  a  trigger  signal  for  the  SSS  routine.  Upon  the  operator's  command, 
the  SSS  routine  is  enabled,  and  the  next  positlve-to-  negative  oroasing 
of  the  filtered  motion  signal  triggers  a  quartz  Crystal  olook/timer. 
Ten  oyolea  of  the  oorreoted  Fy(t)  signal  are  sampled,  at  a  rate  of  too 
samplaa/oyole.  The  second  poaitive-to-negative  crossing  of  the  filtered 
notion  signal  triggers  the  tiner  and  initiates  the  sampling  of  ben 
oyolea  of  the  Fx<t)  signal.  Finally,  the  third  poaitive-to-negative 
oroasing  triggers  the  timer  again,  and  ten  oyolea  of  the  oorreoted  X(t) 
signal  are  sanpled,  Thus,  at  every  test  oondltlon,  1000  data  points  are 
obtained  for  FxU* ) ,Fy(ti ) ,  and  X(ti),  and  the  data  arrays  are  stored  in 
computer  memory. 

Some  Important  points  need  to  be  stressed  oonoerning  this 
foroe/motlon  data  acquisition.  First,  the  bandpass  filter  is  used  only 
to  provide  a  steady  signal  to  trigger  the  tlmer/olook.  Any  modulation 
of  the  motion  signal  due  to  rotor  runout  is  eliminated  by  this  filter, 
aa  long  aa  the  rotational  frequenoy  and  shaking  frequency  do  not 
ooinolde.  Therefore,  the  shaking  frequencies  are  seleoted  bo  avoid 
coincidence  with  running  apeeda.  However,  the  rotor  motion  and  oorreoted 
force  signals  whioh  are  sampled  and  captured  for  ooefriolont 
identification  ore  filtered  only  by  a  low-pass  filter  (500  Hz  cutoff), 
and  the  effects  of  runout  aa  well  as  shaking  motion  a ro  present  In  the 
reoorded  data.  A  second  point  worth  noting  is  that  the  sample  rate  Is 
directly  dependent  on  the  shaking  frequency.  Aa  the  shaking  frequenoy 
is  inoreaaed,  the  sample  rate  (aamples/seoond)  also  increases.  In  order 
to  get  the  desired  100  aamplas/oyole,  shaking  frequencies  must  be  ohosen 
to  correspond  to  disorete  sample  rates  whioh  are  available.  Henoe,  the 
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frequency  at  which  the  rotor  is  shaken’is  oarefully  ohosen  to  provide 
the  desired  sampling  rate  and  a  steady  trigger  signal. 

Moat  of  the  fluid  flow  data  are  used  for  the  input  parameters 
required  by  Nelson's  analysis.  The  upstream  (reservoir)  pressure  and 
temperature,  downstream  (sump)  pressure,  and  the  inlet  circumferential 
velocity  (determined  as  outlined  earlier)  are  provided  directly.  The 
friotion-faotor  values  of  the  rotor  and  stator  are  supplied  in  the 
form  of  Hlrs  coefficients,  which  are  obtained  from  the  pressure 
distribution  data  in  the  manner  described  below. 

Recalling  the  disousalon  of  Hirs'  turbulent  bulk-flow  fluid  model, 
the  model  assumes  that  the  wall  shear  stresses  can  be  written  as  in 
equation  (4).  For  the  gas  seals  discussed  here,  an  adlabatio, 
oompreasible  flow  with  friatlon  analysis  is  employed,  and  the  measured 
pressure  gradient  and  maas  flow  rate  (leakage)  data  are  used  to 
calculate  a  friction  factor  coefficient,  1,  for  eaah  test  condition. 
From  the  1  versus  Rt  and  u>  data,  the  Hirs  coefficients  mr,  nr  of  the 
friction  faotor  formula 

1  -  nr  Ra"P  [  1  +  1  /  4b*]0+»P)/2  ,  b  -  U  /  Ru  (20) 

are  calculated  on  a  least-square  basis.  For  the 
smooth-rotor/smooth-stator  combination,  the  values  are  assumed  to  apply 
for  both  the  rotor  and  stator.  Hence,  for  this  case,  mr-ms  and  nr-ns. 

For  the  amooth-rotor/honeyoomb-stator  combination,  a  combined  A  is 
measured,  which  is  related  to  the  rotor  Ar  and  (honeyoomb)  stator  As  by 

A0  -  (AP  ♦  A0)  t  2  (21 ) 

and  henoe, 

As  -  2A0  -  AP  (22) 


Therefore,  X3  la  determined  from  measured  data  for  X0  and  a 
calculated  value  for  Xr  from  Equation  (21 )  with  experimentally 
determined  values  for  mr  and  nr.  Then,  as  before,  the  Xs  versus  Ra  and  <d 
data  are  used  to  oaloulate  the  Hlrs  coefficients  for  the  honeyoomb 
stator. 

As  stated  previously,  the  Hlrs  coefficients  for  the  seal  rotor  and 
stator  are  required  input  parameters  for  Nelson's  analysis,  as  are  the 
fluid  flow  conditions  up  and  downstream  of  the  seal  and  the  rotational 
speed  of  the  rotor.  The  appropriate  input  parameters  for  each  specif io 
test  case  oan  be  provided  for  Nelson's  analysis  from  statio  test  results 
and  measurements.  In  this  manner,  a  point -by- point  comparison  of 
theoretical  predictions  to  experimental  results  oan  be  made  for  leakage 
through  the  seal,  axial  -  pressure  distribution,  entrance  -  loss 
coefficient,  and  rotordynamio  coefficients. 


TEST  PROCEDURES 


At  *  the  start  of  eaoh  day's  testing,  the  foroe,  pressure,  and 
flowmeter  systems  are  calibrated.  The  total  system,  from  transducer  to 
computer,  is  calibrated  for  eaoh  of  these  variables.  The  foroe  system 
calibration  utilizes  a  system  of  pulleys  and  known  weights  applied  in 
the  X  and  Y- directions.  An  air-operated  dead-weight  pressure  tester  is 
used  for  pressure  system  calibration,  and  flowmeter  system  calibration 
is  aohieved  with  an  internal  precision  quartz  clock  which  simulates  a 
known  flowrate. 

All  of  the  tests  performed  to  date  have  been  made  with  the  rotor 
executing  small  motion  about  a  centered  position.  A  typioal  test  begins 
by  centering  the  seal  rotor  in  the  stator  with  the  Zonlo  hydraulic 
shaker,  starting  airflow  through  the  seal,  setting  the  rotational 
speed  of  the  rotor,  and  then  beginning  the  shaking  motion  or  tne  rotor. 
Data  points  are  taken  at  rotational  speeds  of  200,  500,  and  1000-8000 
opm,  in  1000  opm  increments.  At  eaoh  rotational  speed,  the  Inlet 
pressure  is  varied  and  data  points  are  taken  at  one  unchoked  flow  and 
four  choked  flow  conditions.  For  eaoh  test  oase  (i.e.,  one  particular 
running  speed,  shaking  frequency,  inlet  pressure,  and  prerotation 
condition),  the  measured  leakage,  rotordynamlo  ooef f ioients ,  axial 
pressure  distribution,  and  entranoe  loss  oooffioient  are  determined  and 
reoorded. 

This  test  sequence  is  followed  for  each  of  three  different  shaking 
frequencies,  and  for  three  inlet  3Wirl  directions  (with  rotor  rotation, 
opposite  rotation,  and  no  rotation).  Therefore,  fifty  data  points  are 
taken  per  test  (i.e.  one  shaking  frequenoy  and  inlet  swirl  combination), 
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with  a  total  of  nine  teats  (for*  small  motion  about  a  oentered 
position)  made  per  seal. 
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RESULTS 


The  teat  resulfca  reported  here  were  developed  as  part  of  an 
extended.  Joint  NASA-USAF  funded  reaearoh  program  for  annular  gaa  seal 
studies.  Tests  were  oarried  out  on  a  araooth-rotor/smooth-stator  seal 
and  a  smooth-rotor/honeycomb-stator  seal.  The  dimensions  and  pertinent 
data  for  eaoh  are  included  in  Table  1.  The  honeycomb  stator  insert, 
provided  by  the  Rooketdyne  division  of  Rockwell  International,  is  the 
turbine  interstage  seal  of  the  HPOTP  (High-Pressure  Oxygen  Turbopump)  of 
the  SSME  (Spaoe  Shuttle  Main  Engine).  Fig.  19  illustrates  the  honeyoomb 
configuration. 

The  test  program  had  the  following  objectives: 

1)  Aoquire  leakage,  friotion  faotor,  and  entrance-loss  data  for 
smooth  and  honeyoomb  seals. 

2)  Compare  predictions  from  ourrent  theory  to  test  results. 

3)  Compare  the  stability  performances  of  a  smooth-rotor/ 
honeyoomb-stator  and  a  smooth  oonatani -clearance  annular  seal. 

The  Hirs  coefficients  for  both  seMn  were  determined  In  the  manner 
described  previously.  The  values  of  these  coefficients  are  listed  in 
Table  2.  Relative  roughness  values  based  on  measured  Hirs  ooerrioients 
as  suggested  by  Colebrook  [22]  are  also  included.  Colebrook's 
formulation, 

(4  nsRa“a)-1/2  -  -2  log  ( ( ( e/2C )  /  3.7)  *  (2.51  /  (R^  nsRams  )))  (23) 

is  a  curve-fit  of  experimental  data  obtained  for  fluid  flow  through 
pipes  with  various  wall  roughnesses.  The  appropriate  stator 
coefficients  are  substituted  to  obtain  the  relative  roughness  (e/2C) 
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Table  1.  Teat  seal  specif loationa. 


Smooth-rotor/ 
smooth  stator 

Rotor 

i 

Diameter;  15.136  oh  <5.959  in) 

Material i  304  Stainless  steel 
Surfaoe 

roughness;  0.102  ym  (4  yin) 


Smooth-rotor/ 
honeycomb  stator 

14.453  cm  (5.690  in) 
304  Stainless  steel 

0.127  yin  (5  yin) 


Stator 

Diameter;  15.283  om  (6.017  in) 

Material;  304  Stainless  steel 
Surfaoe 

roughness;  0.140  ym  (5.5  yin) 


14.614  om  (5.754  in) 
6061-T6  Aluminum 

1.575  mm  (0.062  in) 
Comb 


Radial  olearance:  0.7366  mm  (29  mil)  0.8065  mm  (31.6  mil) 

Seal  length;  5.080  00  (2.00  in)  2.540  om  (1.00  in) 


Table  2.  Frlotion-faotor  data. 


Smooth-rotor/ 

smooth-stator 

Smooth-rotor/ 

honeyoomb-stator 

Mira  ooeffioients 

Rotor  ns t 

i 

0.187 

0.187 

Rotor  mat 

-0.333 

-0.333 

Stator  ns: 

0.187 

0.187 

Stator  ms: 

-0.333 

-0.0778 

Relative 
roughness  c/2C 

Rotor: 

1 .44x10“^ 

1  .UUxlO”14 

Stator i 

1  .iiMxlO”*11 

4.93X10-^ 

(Notei  The  relative  roughness  values  shown  are  averages  over  an  axial 
Reynolds  number  range  of  20,000  -  80,000.) 


values.  It  should  be  noted  that  friction  faotor  data  for  honeyoomb 
seals  have  not  been  previously  published.  , 

The  results  provided  here  are  grouped  in  statio  (leakage,  pressure 
distribution,  entrance  loss  coefficient)  and  dynatnio  (rotor dynamic 
coefficient)  sections.  A  one-to-one  comparison  of  the  smooth  and 
honeyoomb  seal  configuration  is  precluded,  however,  due  to  differences 
in  seal  length,  nominal  clearance,  and  inlet  guide  vane  configuration, 

as  illustrated  in  Pigs.  12  and  14. 

\ 

Statio  Results.  Figs. 20  and  21  and  Table  3  show  a  comparison  of  the 
theoretical  and  experimental  leakage  through  the  seal  for  various  fluid 
prerotation  conditions.  The  figures  show  the  leakage  at  various 
proesure  ratios  (reservoir  pressure  /  sump  pressure).  In  the  table,  the 
leakage  has  been  averaged  over  all  speeds  and  pressure  ratios,  and  is 
presented  in  ratio  form  (Theory/Experlment ) .  The  comparison  shows  that 
for  both  the  smooth  and  the  honeyoomb  seal,  leakage  is  underpredioted 
for  the  no n- prerot at ad  case.  Conversely,  for  fluid  prerotation  either 
in  or  opposing  the  direction  of  rotor  rotation,  the  leakage  is 
overpredloted  for  both  seals.  The  maximum  error  is  approximately  7.5%, 
ooourring  for  the  smooth  seal  with  prerotation  in  the  dlreotion  of  rotor 
rotation.  Average  error  for  the  smooth  seal  is  1.7#,  and  for  the 
honeyoomb  seal  is  1#. 

The  pressure  gradient  plots  are  included  in  Figs, 22-28.  Fig.  22 
illustrates  the  negligible  effoot  of  running  speed  on  the  pressure 
distribution  in  the  seal.  This  plot  has  ten  curves  (corresponding  to 
the  ton  rotational  speed  increments)  plotted.  This  accounts  for  the 
heavy  lines  whioh  appear  in  some  oases.  This  particular  plot  is  of 
the  experimental  data  for  the  non-prerotated  smooth  seal  case,  however, 
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Table  3.  Theory  versus  experiment  leakage  comparison. 
(Thaory/Experiment ) 


Fluid  prerotation 
dlreotion 

Smooth-rotor/ 

smooth-stator 

Smooth-rotor/ 

honeyoomb-stator 

With  hotor  rotation 

1,075 

1.047 

(standard  deviation) 

0.012 

0.008 

No  prerotatlon 

0.9684 

0,9712 

(standard  deviation) 

0.012 

0.005 

Opposite  rotor 
rotation 

1.007 

1.013 

(standard  deviation) 


0.012 


0.004 


pressure  gradient  in  saooth  seal  with  no  prerotation,  unchoked  flow, 
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none  of  the  pressure  plots  show  any  appraolabls  variation  due  to  running 

spaed. 

Figs.  23-28  show  the  theoretloal  and  experimental  pressure  data  for 
eaoh  of  the  seals  wider  various  prerotation  conditions.  Due  to  the 
absenoe  of  running  speed  dependence,  only  one  speed  is  plotted  for  each 
Inlet  pressure  condition.  The  numbers  on  the  plotted  lines  refer  to  the 
Inlet  pressure,  where  1  corresponds  to  approximately  166.9  kPa  (12.H 
psig),  and  2  through  5  correspond  to  308,2  kPi  (30  pslg),  446,1  kPa  (50 
psig),  584.0  kPa  (70  psig),  and  721.9  kPa  (90  .psig),  respectively,  The 
lowest  pressure  corresponds  to  unohoked  flow  through  the  seal,  while  the 
others  are  choked.  The  shapes  of  th<>  pressure-gradient  plots  show 
fairly  good  correspondence  between  theory  and  experiment.  This  to  be 
expected,  however,  slnoe  the  Hire'  coefficients  used  in  the  analysis 
aome  dlreatly  from  the  measured  pressures.  The  best  agreement  occurs 
for  the  non-prerotated  flow  in  both  the  emooth  and  honeyoomb  seals.  For 
prerotated  flow  in  either  direction,  the  theoretloal  gradient  Is  shirted 
up  slightly  for  both  seal  configurations.  This  upward  shift  is  partly 
due  to  a  total  preaaurs  oorreotlon  that  is  made.  When  the  now  is 
prerotated  by  the  guide  vanes,  it  is  aooelsrated  an  well  as  turned,  and 
the  measured  statlo  pressure  at  the  vane  exit  decreases.  This  explains 
why  the  experiments!  plots  show  lower  seal  entrance  pressures  Tor  either 
prerotated  case  than  for  tha  non-prerotated  oeee.  Nelson's  analysis, 
however,  assumes  thst  the  supply  pressure  upstream  of  the  seal  is  the 
total  pressure.  Hanot,  the  axial  oomponent  of  the  fluid  veloolty  as  it 
leaves  the  guldo  vanes  is  used  to  oalculate  an  effeotlve  total  pressure, 
whioh  is  higher  then  the  measured  static  pressure.  This  corrected 
pressure  is  then  input  as  the  reservoir  pressure  to  Nelson's  analysis. 


sure  gradient  in  honeycomb  seal  with  prerotation  in  the  direction  of  rotor  rotation. 
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The  entranoe  loss  coefficient,  7,  Also  nay  havs  sons  bearing  on  the 
upward  shift  seen  in  some  of  the  pressure  gradient  plots.  Plots  of 
(7*1)  versus  axial  Reynolds  number  are  inoluded  In  Pigs.  29-34.  Ten 
experimental  ourves,  corresponding  to  running  speed  are  plotted  on  eaoh. 
Recalling  Deissler'a  curve  fit  employed  by  Nelson  and  plotted  in  Fig.3> 
experimental  results  show  loss  terms  (7*1 )  significantly  higher  than 
those  predloted  for  both  prerotated  honeyoomb  seal  oases.  Agreement 
between  theory  and  experiment  is  fairly  good  for  the  non- prerotated 
oases  for  both  seals.  For  the  smooth  seal  with  prerotation  in  the 
dlreotion  of  rotor  rotation,  the  loss  ooeffioient  is  overpredioted,  with 
the  experimental  results  indioating  a  negative  7. 

Dynamlo  Results.  Dynamlo  teats  were  performed  at  shaking  frequencies  or 
58.8,  74.6,  and  124.6  Hz.  As  was  dlsouased  in  the  Data  Acquisition 
section  of  this  report,  these  frequencies  were  ohosen  to  provide  the 
desired  sample  rate  and  a  steady  trigger  signal.  The  dynamlo 
coefficients  obtained  at  the  two  lower  frequencies  are  essentially  the 
same.  At  the  124.6  Hz  shaking  frequenoy,  however,  correspondence  of  the 
data  to  that  obtained  at  the  lower  frequencies  is  unsatisfactory.  In 
seeking  to  explain  the  disorepanoy,  tests  were  run  to  determine  the 
relative  transfer  funotion  of  the  teat  apparatus.  The  plots  in  Fig.  35 
show  the  results  of  these  tests,  and  indicate  a  reaonanoe  of 

the  apparatus  occurring  at  approximately  25  Hz  (the  drop  in  phase 

difference  at  approximately  45  Hz  corresponds  to  a  resonance  of  the 
shaker  support  structure).  As  the  shaking  frequenoy  is  increased  above 
this,  the  input  force  levels  required  to  aohieve  a  given  motion 
amplitude  increase  rapidly.  At  the  124.6  Hz  shaking  frequency, 

attainable  motion  amplitude  is  about  50%  of  that  achieved  at  the 


Entrance-loss  for  honeycomb  seal,  no  prerotatlon 
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58,8  and  74.6  Hu  frequencies.  Therafore,  ona  possible  explanation  for 
the  poor  agreement  between  the  results  is  that  as  motion  amplitude 
decreases,  so  does  the  foroe  measured  by  the  load  oells,  and  the 
measurement  system  resolution  suffers. 

The  plots  of  the  rotordynamio  coefficients  are  found  in  Pigs. 
36-43.  These  plots  include  both  the  theoretical  and  experimental  data. 
The  coefficients  are  plotted  versus  the  reservoir  t  sump  pressure  ratio, 
and  the  solid  lines  correspond  to  the  theoretioal  data.  The 
experimental  results  are  indicated  by  symbols.  The  location  of  the 
symbols  represents  the  average  value  of  the  ooeffiolent  (averaged 
over  all  of  the  running  speeds)  at  each  particular  inlet  pressure 
condition,  and  the  vertical  lines  drawn  through  the  symbols  signify  the 
range  over  whioh  they  varied  through  the  speed  range,  The  test  results 
plotted  here  were  obtained  by  shaking  the  rotor  with  an  amplitude  of 
seven  mils  at  74,6  Hz. 

Dynamio  Results  -  Smooth  Seal.  For  the  smooth  seal,  direct  stiffness 
(Fig. 36)  is  overpredioted  for  the  non-prerotated  condition,  and 
under predicted  for  both  prerotated  conditions.  Best  agreement  is  seen 
in  the  oase  for  prerotation  opposite  rotor  rotation,  and  the  trend  of 
inoreaaing  stiffness  with  increasing  pressure  ratio  compares  favorably. 

In  the  oross-aoupled  stiffness  comparison  for  the  smooth  seal 
(Fig.  37),  theory  overpredicts  for  both  prerotation  conditions,  and 
underprediots  for  the  straight  flow  case.  In  every  instanoe,  agreement 
becomes  worse  with  increasing  pressure  ratio.  The  non-prerotated  case 
shows  a  divergence  both  in  magnitude  and  sign.  It  should  be  noted, 
however,  that  the  magnitudes  for  this  case  are  significantly  smaller 
than  for  either  of  the  prerotated  oases. 
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Fig.  37  Cross-coupled  stiffness  of  ssooth  seal 
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Agreement  between  theory  and  experiment  for  direct  damping  (Fig, 38) 
is  the  moat  favorable  of  all  the  dynamic  coefficients.  Direct  damping 
for  the  smooth  seal  shows  an  increase  for  increasing  pressure  ratio 
aoross  the  seal,  with  theory  over predicting  slightly  for  the  case  of 
prerotation  in  the  direction  of  rotor  rotation.  For  the  other 
prerotation  conditions,  the  direot  damping  is  slightly  overpredicted  at 
higher  pressure  ratios. 

Cross-coupled  damping  (Fig, 39)  for  the  smooth  seal  generally  shows 
agreement  in  the  trends  for  the  theoretical  and  experimental  results. 
For  prerotation  in  and  opposing  the  direotion  of  rotor  rotation,  the 
theory  underprediots  oross-coupled  damping  magnitude  by  approximately 
509,  but  shows  a  sign  consistent  with  the  test  data.  For  the 
non-prerotated  case,  the  theory  predicts  coefficients  so  small  as  to  be 
considered  negligible.  This  is  not  inconsistent  with  the  test  results, 
however,  as  the  magnitudes  for  this  oase  are  significantly  smaller  than 
for  either  prerotated  case. 

Dynamic  Results  -  Honeycomb  Seal.  The  honeycomb  seal  data,  in 
general,  shows  the  same  correspondence  between  theory  and  experiment  as 
the  smooth  seal.  A  notable  exception,  however,  is  in  the  direct 
stiffness  coefficient  (Fig. 40).  For  both  prerotated  cases,  the 
predicted  stiffness  decreases  with  increasing  pressure  ratio,  while  the 
measured  stiffnesses  Increase.  This  same  predicted  decreasing  trend  is 
shown  for  the  non-prerotated  case  at  the  lower  pressure  ratios, 

In  the  oross-coupled  -;iffness  comparison  (Fig.4l),  the  theory 
underprediots  the  magnitudes,  but  correctly  predicts  the  signs  of  the 
coefficients.  For  the  non-prerotated  oase,  the  predicted  stiffnesses 
are  essentially  zero.  The  relative  magnitudes  of  the  experimental 


results  for  this  same  case  In  comparison  to  either  prerotated  case  are 
also  quite  small,  however.  For  prerotation  in  the  direction  of  rotor 
rotation,  theoretical  cross-ooupled  stiffnesses  are  approximately  25 1> 
less  than  experimental  ones.  For  oounter  prerotation,  theory 
underprediots  by  about  50* . 

With  the  exception  of  the  non- prerotated  case,  agreement  between 
theory  and  experiment  is  fairly  good  for  the  direot  damping  coefficients 
(Fig. 42)  of  the  honeyoomb  seal.  In  the  non-prerotated  oase,  theory 

underprediots  the  ooeffioients  by  approximately  46*.  The  prerotated 
oases  show  agreement  to  within  approximately  10*. 

Theoretical  results  for  the  oross-coupled  damping  ooeffioients 
(Fig. 43)  of  the  honeyoomb  seal  are  small  enough  to  be  considered 

negligible.  In  every  oase,  the  theory  underprediots  the  ooeffioients  by 
a  wide  margin.  However,  the  trend  of  increasing  magnitude  with 

increasing  pressure  ratio,  as  well  as  the  signs  of  the  ooeffioients, 

agree. 

One  method  in  whioh  the  dynamic  ooeffioients  of  the  smooth  and 
honeyoomb  seals  oan  be  direotly  compared  is  through  their  respective 
non-dimensional  whirl  frequency  ratios.  Whirl  frequency  ratio  is  defined 
Whirl  frequency  ratio  -  k  /  CJl  , 

where  ft  is  the  shaking  frequency.  This  ratio  is  a  measure  of  the 
destabilizing  influence  of  the  cross-ooupled  stiffness  with  respect  to 
the  stabilizing  influence  of  direot  damping.  Plots  of  whirl 
frequency  ratio  versus  running  speed  are  included  in  Fig. 44.  The  smooth 
seal  plot  shows  a  small,  positive  whirl  frequency  ratio  over  most  of 
the  running  speed  range.  The  honeycomb  seal  plot,  however,  shows  a 
negative  whirl  frequency  ratio.  The  negative  sign  arises  due  to  a 
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negative  oross-coupled  stiffness.  This  negative  k  exerts  a  stabilizing 
influence,  resulting  in  a  force  which  acts  in  the  same  direction  as  the 
direct  damping  force. 
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CONCLUSIONS 


A  seal-test  facility  has  baan  developed  for  the  study  of  various 
types  of  gas  seals.  A  method  of  determining  rotordynamio  ooeffioients 
has  been  established,  and  consistent,  repeatable  results  have  been 
obtained.  After  some  initial  failures  in  the  .test  apparatus, 
reliability  has  been  satisfactory,  and  a  oomplete  set  of  experimental 
results  can  be  acquired  in  a  matter  of  days. 

The  experimental  and  theoretical  results  of  the  preceding  section 
support  the  following  oonolusionst 

(a)  Theoretical  results  for  leakage  are  consistent  with  test 
results.  Slightly  higher  leakage  ocours  for  oases  where  the  flow  is  not 
prerotated.  Agreement  between  theory  and  experiment  is  satisfactory, 
with  a  maximum  error  of  approximately  7.5%. 

(b)  Experimental  and  theoretical  results  for  the  pressure 
distributions  and  entranoe-loss  coefficients  are  relatively  Insensitive 
to  running  speed  for  the  ranges  (0-8500  opm)  and  seals  tested  to  date. 

(c)  The  entrance-loss  relationship  (Deissler  [13])  employed  by 
Nelson  is  inconsistent  with  the  test  results  for  oases  where  the  fluid 
is  prerotated.  For  all  auoh  cases  exoept  one,  the  entrance-loss 
coefficient  is  underpredicted. 

(d)  In  the  test  results  for  the  honeycomb  seal,  the  steep  entrance 
pressure-loss  seems  to  extend  partially  inside  the  seal.  Also,  the 
measured  pressure  at  the  exit  of  the  seal  generally  equals  the  back 
pressure,  rather  than  being  greater,  as  is  predicted  by  compressible 
flow  theory  for  choked  flow.  These  same  phenomena  do  not  occur  for  the 
smooth  seal,  implying  that  perhaps  the  effective  length  of  the  honeycomb 


seal  is  less  than  its  actual  physical  length. 

(e)  Test  results  for  the  direct  stiffness  of  both  the  smooth  and 
honeyoomb  seals  show  muoh  greater  sensitivity  to  fluid  prerotatioh  than 
predicted  by  theory.  Prerotation  of  the  fluid  (in  either  direotion) 
results  in  measured  direct  stiffnesses  whioh  are  significantly  larger 
than  for  no  prerotation.  Theory  prediots  the  direot  stiffness  to  be 
relatively  insensitive  to  fluid  prerotation. 

(f)  Theoretical  predictions  of  the  Influence  of  fluid  prerotation 
on  oroso-ooupled  stiffness  and  damping  are  consistent  with  the  test 
results.  In  general,  theory  underprediots  the  magnitudes  of  these 
oross-coupled  coefficients,  while  oorrectly  predicting  their  trends 
with  respeot  to  prerotation. 

(g)  Agreement  between  theory  and  test  results  for  the  direct 
damping  coefficients  is  favorable.  For  both  the  smooth  and  honeyoomb 
seal,  dlreat  damping  is  largest  for  no  fluid  prerotation. 

(h)  Over  the  speed  range  tested,  none  of  the  rotordynamic 
coefficients  show  appreciable  sensitivity  to  the  rotational  speed  of  the 
rotor.  This  may  be  due  to  the  lack  of  development  of  significant  shear 
foroes  in  the  seal.  It  appears  that  running  speeds  above  those  attained 
to  date  may  be  necessary  to  produce  significant  shear  force  afreets. 

(i)  For  the  non-prerotated  case,  the  smooth  seal  has  a  positive 
aross-ooupled  stiffness,  while  k  for  the  honeyoomb  seal  is  negative. 
This  negative  cross- coupled  stiffness,  and  hence  negative  whirl  ratio, 
indicates  that  the  stability  performance  of  the  honeyoomb  seal  is 
more  favorable  than  that  of  the  smooth  seal. 
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ABSTRACT 


A  nav  varalon  of  a  *virl-flov  flnlea  dlffacanoa 
coaputar  program  waa  davalopad  lsoorpotatlni  tba  recent 
quadratic,  upwind-weighted  differencing  tcheae.  It  waa 
aaployad  In  pradlctln|  tba  lncoapraaalbla  flow  within 
aach  of  two  ilxpla  labyrinth  aaal  cavity  eonflgura- 
tlona.  Cloaa  agreement  batwaan  pvadietlona  and  pre¬ 
vious  comap  ending  aaal  mass ur  amenta  la  abown  for 
verification  puxpoaaa.  Tba  dlatrlbutlon  of  aach 
quantity  revealing  tba  complete  flow  ebaraotar  for  a 
labyrinth  aaal  cavity  la  praa anted  and  dlaouaaad  bora 
for  tba  flrat  time.  Tba  affact  of  aaal  cavity  geometry 
on  overall  praaaura  drop  la  abown  to  be  significant. 

In  addition,  thla  lnvaatlgetlon  conatltutaa  a 
timely  evaluation  of  tba  QUICK  convactiva  differencing 
■ebana  for  raallatlc,  highly  non-llnaar  turbulent  flow 
applleatlona.  No  algnlflcant  atablllty  dlfflcultlaa 
wara  encounter  ad.  Ho  mover,  the  r  faulting  increased 
accuracy  allow*  an  accurate  aolutlon  uaiag  conaldarably 
fewer  grid  point*  than  for  tho  Hybrid  upvlnd/caneral 
scheme,  Computer  execution  tin*  con— mad  aaploylng 
QUICK  waa  only  44X  of  that  for  tba  Hybrid  arhana. 
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Axial  length  of  aaal  cavity 

Grid  Pdclet  nuahar 
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Axial  taynolda  number 
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Bulk  awlrl  velocity 
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c 

Turbulanea  dlaaipatlon  rata 

* 

Ganarallaad  dependant  variable 

p 

Til— -naan  density 

V 

r in— atlc  viscosity 

n 

Shaft  epeod 

Subscript* 

n.a.a.w 

North,  south,  a— t,  wait  face*  of  a  con¬ 
trol  voluna 

H.S.I.W 

North,  South,  East,  Wait  neighboring 
grid  points 

OmODOCTlON 
jifllhg  yputtd 

Labyrinth  aaal*  play  a  vital  tola  In  gna  and  ataau 
turblnca,  aoapraaaora,  and  hlgh-capaclty  fuel  puapa. 

In  fuel  puapa  tba  aeellag  objective  la  generally  to 
nlAlnlaa  tba  leakage  flew  around  the  lapallara,  wbaraaa 
for  tucblaaa  It  la  oftan  to  naxntaln  a  minimal  leakage 
flow  for  cooling  hot  part*.  The  baale  idea  in  tbs 
design  of  a  labyrinth  aaal  la  to  ptovlda  a  highly  dla- 
ilpatlve  flow  path  batwnau  high  and  low  praaaura  re¬ 
gion.  Thla  la  acconpllahad  by  naan*  of  a  earies  of 
restrictor  teeth  and  dissipation  cavltiaa.  Tba  flow 
paaaaga  through  a  labyrinth  la  llluatratad  In  fig,  1 
for  a  a lap la  atralght- through  aaal  configuration.  Each 
tooth  convert*  a  portion  of  the  available  praaaura  head 
into  naan  flow  kinetic  energy,  soma  of  which  la  dls- 
elpated  within  the  cavity  1— dlataly  dowuatraaa. 

Hunt roue  lnvaatlgatora  have  mxperimantelly  atudlad 
several  groaa  overall  flow  charactarlatles  of  labyrinth 
seals  la  attainting  to  obtain  empirically-baaed  r ala- 
dona  for  estimating  tba  laakaga  rata,  A  few  have 
aaployad  lncoapraaalbla  fluids  [1-A]  while  others  have 
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utilised  coxpreeeible  onee  [5-10],  of tan  at  ineompree- 
aibla  flak  condition*.  Tha  majority  of  thaaa  investi¬ 
gation*  hi.va  performed  parametric  studies  concerning 
how  pressure  drop,  aaal  design,  claarancaa,  rotational 
tpaad,  ate.  affaet  overall  performance.  Tba  raault* 
of  euch  studies  can  ba  uaad  with  a  fair  amount  of 
accuracy  for  predicting  the  leakage  rata  through  tha 
particular  labyrinth  aaal  that  we  atudied.  However, 
whan  a  labyrinth  aaal  of  a  eonaidarably  different 
dailgn  la  conaldared,  uae  of  thaaa  smplrleally-baasd 
relatione  antalla  aubatantlal  error. 

Accurata  nuaarical  simulation  of  labyrinth  aaal 
flowflalda  smploying  the  governing  differential 
aquation*  la  relatively  complicated  and  expansive. 

A  recant  computerised  literature  aaarch  conducted  by 
tha  praaast  authora  Indicated  that  oaly  Itoff  [11]  end 
Rhode  at  al.  [12]  have  reported  aueh  an  Inveatigatlon. 
Scoff's  eolation  corresponded  to  a  very  large  aeala 
experimental  water  facility  eonaiatlag  of  a  etralght- 
through  aerlea  of  generic  eevitiae  aa  illuatratad  in 
rig.  1.  Ha  eoapared  a  aingla  radial  profile  of  pre¬ 
dicted  naan  evlrl  velocity,  r.rn.e.  twirl  velocity  -nd 
turbulence  dlaalpatlon  rata  with  hie  correaponding 
maaauramanta  for  an  axial  atation  Midway  between 
adjacent  taath.  A  sample  of  theaa  maaeur amenta  waa 
utlliaed  la  thla  Inveatigatlon  for  experimental  veri¬ 
fication  purpoaaa  (thaaa  are  apparently  tha  only 
detailed  measurements  available).  Rhode  at  al.  pra- 
aantad  aaveral  aolutiona,  involving  different  grida 
and  convection  differencing  achcnna,  dearly  illua- 
tratlng  fdae  dlffualon  nuaarical  error  effacta  on  the 
aolutlon  for  compraaalble  flow  through  a  cavity  of 
unique  configuration. 


Stator  Housing 


jianonoflofL 


Rotating  Surface 
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Fig.  1  Generic  labyrinth  aaal  configuration 


Felee  Dlffualon 

Until  recently  tha  wort  effective  method  for 
avoiding  laetabillty  in  numerical  aolutiona  of  con¬ 
vection-dominated  flowa  hma  been  tha  uae  of  upwind 
differencing  for  convection  tana.  However,  thla 
flret-order  accurata  achewe  Introduce*  felee  dlffualon 
[13].  Thla  la  a  eacond -order  truncation  numerical 
error  which  can,  under  certain  condition* ,  yield 
eerioua  conaequencaa  for  aolutlon  accuracy.  Thla  dif¬ 
ficulty  occur*  where  convection  dominate*  (l.a.  where 
the  grid  Pdclet  number  jle|  ■  |v|a/1>  exceed*  2.0)  in 
the  preaenca  of  both  stream) Ine-to-grld  ikewnaaa  and 
diffualve  traaaport  normal  to  the  flow  direction. 

Pale*  dlffualon  raault*  In  aa  overly  diffualve 
aolutlon.  la  fact,  acme  aolutiona  have  been  obtained 
In  which  fdae  dlffualon  haa  entirely  obacured  phye- 


ical  dlffualon  [14-15].  It  la  the  turbdent  recir¬ 
culating  type  of  flow,  euch  aa  that  of  the  preeant 
problem,  which  la  particularly  auacaptible  to  thla 
phenomenon.  Grid  refinement  eufficlent  to  allow  can- 
trd  differencing  (l.e.  whan  | Fa |  i  2.0)  for  convec¬ 
tion  tnraa  would  clrcuxvant  thla  error.  However,  coat 
and  atoraga  dlffieultlaa  would  bacona  aevara  for  moat 
two-dimanolond  problem*  and  prohibitive  for  three- 
dlmenalond  ones. 

Fortunately,  two  improved  convective  differencing 
eppvoecbea  which  miaiaiae  fdae  dlffualon  have  vacant¬ 
ly  baan  propoaad.  Tha  "akaw  upwind  differencing" 
aehaae  of  Kalthby  [16]  raducea  akawnaaa  error*  by 
aaaactlally  aligning  tha  upwind  diffarennlag  acheme 
for  convection  tarna  along  local  atraamlinaa.  Tha 
recant  QUICK  (quadratic  Upstream  Interpolation  for 
Convactiva  Kinematic* )  scheme  of  Leonard  [17]  eliai- 
aataa  fdae  dlffualon  by  incorporating  a  three-point 
upwind-ehlfted  interpolation  formula.  lecant  in¬ 
vestigation*  by  Latch* Iner  [IS]  and  Leech* inar  and 
Rod!  [19]  dlacuat  raault*  using  both  QUICK  and  skew- 
upwind  dlffaranelag,  Indicating  dlatlnot  advantage* 
over  tha  hybrid  upwind/eantral  dif far anting  achama. 

Although  aavard  applluatlon  atudlaa  [18-24] 
have  Incorporated,  one  or  both  of  thaaa  two  recant 
k shame*,  moat  of  tha  flow  problama  conaldared  thua  far 
are  relatively  easy  to  eolva.  Hence  tha  atablllty  and 
ovarahoet  tandanoiaa  of  thaam  two  achomaa  remain  to  ba 
luffleiantly  challangad.  Spmciflcally,  to  tha  knowl¬ 
edge  of  thaia  author*,  only  Laaehalaar  and  Rodl  [19] 
and  Han  at  d.  [22]  have  formally  documented  predic¬ 
tion*  ualag  QUICK  for  highly  non-llnaar,  turbulent 
recirculating  flown.  Sinilsrly  only  Mllitanr  at  d. 
[24],  Laaehalaar  and  Rodl  [19],  and  Baaodakar  at  d. 
[23]  have  reported  using  Kalthby' ■  mathod  for  auch 
caaoa.  Tha  praaact  problaa  la  porhapa  tha  moat 
challenging  for  atablllty  considerations  baaauaa  of 
the  vary  high  la  and  the  enhanced  praaeura-vdoolty 
noa-liaaarlty  due  to  tha  praaaaee  of  awlrl  momentum. 

NUMERICAL  PROCKDOU 

Oenaral  Remark a 

The  numerical  procedure  amployad  in  thla  In¬ 
vestigation  ia  baaed  on  that  of  tha  T1AC1-2I  computer 
program  of  Goaman  and  Pun  [25],  Tha  recant  QUICK 
com  active  differ  too  lag  achama  of  Lao  nerd,  which  ia 
described  below,  wee  Implemented  ee  e  new  version  of 
the  code.  Other  ( naturae  Include  swirl  momentum  and 
ita  Influence  on  wall  function*  an  wall  as  tha 
simulation  of  oloplng  wall  boundaries .  Thar*  are 
aevard  method*  of  numerically  simulating  tha  sloping 
cavity  wall  ahem  la  Pig.  5(b).  In  the  interest  of 
conceptual  simplicity,  tha  atairstap  alnulatlon  ap¬ 
proach  Illuatratad  in  rig.  2  waa  followed.  Rhode  et 
d.  [26]  found  that,  if  carefully  implemented,  the 
error  arising  from  this  approximation  ia  small  in 
comparison  to  that  of  the  turbulence  model. 

The  governing  equations  for  thla  arl symmetric 
flow  consist  of  conservation  of  maaa  end  momentum 
(with  x,  r,  •  velocity  component*  u,  v,  w),  as  mall  as 
transport  equations  for  turbulanea  kinetic  energy  k 
end  turbulence  dlselpatlon  rate  c.  They  ney  be 
expressed  In  tba  general  form 

-  l|y  (pm  ♦)  ♦  It  <P"  -  fe  <*%  - 

fc<*r4  -  *♦  (1) 


Fl|.  2  Computational  domain  for  a  labyrinth  aaal 
cavity  with  a  coaraa  grid 


whare  d  vapcaaanta  any  of  tha  dapaadant  var labial,  and 
tha  aquations  dlffar  primarily  in  chair  final  aeurea 
tana  S4  [27].  It  la  tha  standard  varalon  of  tha  two- 
aquatlon  k-c  turbulanea  modal  [28] ,  axtandad  to  lnoluda 
swirl  mean  turn,  which  boa  baas  aaployad  thua  far.  Tha 
corresponding  finite  dlffar anti  aquation*  ara  aclvad 
via  tha  SIMM  ml* implicit,  llna-by-lina  aathod 
employing  tha  trldlagonal  aatrla  algorithm.  Tha 
aolutlan  aathod  advancaa.  procaaalng  ou  aoluaa  of 
eontrol  voluaaa  at  a  tlaa,  alvaya  ava aping  in  tha 
downstream  dlractlon. 

Boundary  valuaa  at  tha  cavity  Inlat  ara  naturally 
vary  important,  but  unfertunataly  vara  unknown,  aa 
nacaaaary  quantltatlva  flow  aaaauxamanta  wara  unavail* 
abla  avan  at  a  alngla  oparatlng  condition.  Dua  to  thin 
lash  of  Inlat  boundary  data,  for  aanh  numerical  ltara*  - 
tlon  tha  Inlat  valuaa  of  aach  variable  (aaoapt  proa* 
aura)  wara  aaaumad  equal  to  tha  lataat  aai-raapondlag 
outlat  valuaa.  Thla  practice  lapUaa  tha  praaaooa  of 
a  aarlea  of  geometrically  identical  oavltlaa  and  a 
atraamwlao  periodic  flow  occurring  downatraaa  of  tho 
flrat  two  or  thrao  oavltlaa.  Tha  boundary  treatment 
utlllaed  along  eash  wall  aaployad  the  wall  function 
formulae  for  twirling  flow*  [27].  A  apaclal  technique 
of  Implementing  thaaa  formulae  waa  utlllaed  for  con¬ 
trol  volume*  along  tha  rotor  In  order  to  proparly 
account  for  tha  circumferential  velocity  of  tha  aurfaca. 


The  computer  cod*  originally  mabodlad  the  lybrld 
upwind /central  differencing  achama  for  convection  tarna. 
Thla  formulation  ylalda  tha  upwind  dlffaraoea  craatmant 
whan  tha  grid  Moist  number  nagnitud*  In  n  given 
dlraotlon  ancaada  2.0|  otherwise  central  differencing 
la  uaad.  Tha  falaa  dlffualon  difficulty  with  tha 
firat-ovdar  aocurat*  upwind  portion  of  tha  lybrld 
achama  motivated  tha  uaa  of  QUICK.  Tha  QUICK  ochama 
ylalda  a  thlrd-ordrr  accurate  differ east  aquation  which 
la  free  of  falaa  diffusion,  a  second  order  truncation 
error.  Lalthby' •  aksw  upwind  differ anting  achama  would 
have  probably  dispensed  with  falaa  dlffualon  juat  aa 
affectively. 


Following  tha  approach  of  Goaman  and  Fun,  tha 
flnlta  difference  aquation  for  aach  4  can  ba  obtained 
by  first  analytically  integrating  Eq.  1  over  tha  volume 
of  tha  appropriate  micro-control  volume  using  tha  Cause 
Divergence  thaorma.  Tha  result  may  ba  expressed  in 
terms  of  tho  control  volume  fata  valuaa.  Tha  eon- 
vactioo  and  diffusion  terms  become  aurfaca  Integrals 
of  tha  convaotlv*  and  diffusive  fluxes  while  tha  aourca 
tarn  la  linear  la  ad.  This  results  In  [23] 
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where  subscripts  n,  a.  a  and  w  retar  to  north,  south, 
seat  and  waat  control  volume  fates.  There  at*  various 
ways  of  evaluating  the  taca  valuaa  of  d  1*  tha  convic¬ 
tion  terms  using  lataat  available  valuaa  at  neighboring 
nodes.  Each  aathod  of  evaluation  conatltutaa  a  dis¬ 
tinct  conviction  differencing  b chime.  Tha  obj active 
of  any  method  la  to  accurately  evaluate  tha  oonvaeted 
d -value  at  each  face  et  tha  control  voIum  while 
avoiding  rhe  Instability  associated  with  high  convec¬ 
tion  rataa. 

Observe  from  tha  illustration  in  rig.  3  that  tha 
upwind  achama  la  based  on  the  uniform  distribution 


(a)  Upwind  differencing  concept 


(b)  qCXCK  differencing  concept 


Fig.  3  Illustration  of  (a)  Upwind  sad  <b)  QUICK 

dlffar seeing  concept*  for  evaluating  d  at  a 
waat  control  volume  face  vith  flow  from  loft 
to  right 


.assumption  In  evslusting  4  at  any  faca.  Also  shown  la 
that  QUICK  la  based  on  upatraam-ahifts d  quadratic  in¬ 
terpolation  for  aaoh  control  volume  faca.  The  wait 
faca  QUICK  interpolation  expressions  using  a  uni fora 
grid  for  simplicity  ara 

♦„  -  K  *  V  -  J<*P  -  2*w  +  V 

if  uv  la  positive  and 

♦w  "  I(*p  +  V  *  I<*x  “  *♦?  +  V 

If  u  la  nagstlva.  The  first  tax*  In  thasa  expres¬ 
sion!  rapraaanta  the  cantarad  dlffaranca  foraula,  and 
tha  aacond  la  tha  crucial  stabilising  upatraan  weighted 
curvature  contribution.  Quadratic  interpolation  was 
incorporated  only  for  tha  cosvaotad  variable  *  wharaaa 
cantarad  dlffaranclng  was  applied  for  the  convncting 
velocities. 

This  acheaa  has  bean  carefully  implemented  In  tha 
convaetlva  taxma  of  only  tha  eonentun  aquations. 
Laachalnar  and  lodl  [19]  have  supported  this  in  ex¬ 
plaining  that  solutiona  for  k  and  e  ware  found  to  be 
unaffected  by  tha  oonvaotlvn  dlffaranclng  eeheaa. 
aapaclally  In  tha  shear  layer  bordering  a  recircula¬ 
tion  tone.  They  attribute  this  to  tha  feet  that  tha 
source  tana  of  thaaa  aquations  dominate  convective 
transport  vithln  tha  influential  shear  layer  auoh  as 
that  found  In  the  present  aeal  problem. 

RI3ULTS  AMO  DI8CUSSI0H 

gxparimantal  Variflostion 

A  computarlaad  literature  search  irdieatad  that 
only  tha  paper  by  Stoff  [11]  givaa  detailed  quantita¬ 
tive  velocity  measurement!  within  a  labyrinth  seal 
cavity.  Only  the  circumferential  velocity  component 
la  presented!  it  was  measured  with  a  laser  dopplar 
anemometer  at  a  ling la  axial  station,  loostsd  midway 
between  adjacent  teeth.  These  measurements  wars 
employed  for  preliminary  verification  testing  of  tha 
current  computer  code. 

Stoff  utilised  a  large-scale  test  section  modal 
of  a  straight-through  anal  with  raotamgular  cavities 
through  which  water  flowed.  It  is  assumed  that  tha 
flow  conditions  ara  characterised  primarily  by  tha 
axial  Reynolds  number  Imx  ■  2Ce/v  *  3. A  z  10*  and 
tha  Taylor  nuaber  Ta  I  (Wd/vi) (d/r,h)’  %  1.1  s  10*. 
Thaaa  ara  baaed  ant  tooth-clearance  bulk  axial 
velocity  U.  tooth  clearance  c,  Vincis* tic  vieooelty  u, 
cavity  bulk  swirl  velocity  H,  shaft  radius  rah«  cad 
radial  distance  from  cavity  baas  to  stator  wall  d. 

Figure  4  shows  that  predictions  using  tha  com¬ 
puter  ooda  described  above  ara  in  good  agreement  (6* 
discrepancy)  with  Itoff'a  corresponding  measurements. 
Further  k  Staff.'  t  computation  employed  tha  asm*  numer¬ 
ical  approach  and  Hybrid  differencing  scheme  discussed 
earlier,  and  bis  predictions  are  essentially  identical 
to  thoea  ahown  hare. 

9FfIfftfri.Sffi41t4f.ai 

Tha  operating  conditions  of  tha  two  seals  inves¬ 
tigated  herein  ara  identical.  Liquid  hydrogan  at  42*1 
enters  a  cavity  with  tooth-clearance  bulk  velocity 
U  -  338  a/e.  The  primary  dlnaaaionleae  flow  ssra- 
astari  (defined  earlier)  arat  lag  ■  1.0  x  10*  and 
Ta  ■»  1.2  x  10s.  Other  conditions  are)  abaft  apaad 
(1  •  33.410  rpa.  cavity  axial  length  L  ■  1.113  m. 
stator  wall  radlua  1  ■  42.89  am.  tooth  radial 


Fig.  4  Compariton  of  present  predictions  (————■) 
of  dimanaionlaas  twirl  velocity  with  cor¬ 
responding  IDA  measurements  (m,  (  H) 
of  Stoff  [11 ]  at  x/L  ■  0.3 

clearance  c  -  0.216  ms.  radial  distance  from  cavity 
beta  to  tutor  wall  d  -  1.103  mi.  and  klnamatlc  viscos¬ 
ity  y  •  1.462  x  10"’  m*/e. 

A  grid  dependence  study  was  conducted  and  is  dis¬ 
cussed  in  detail  below.  Convergence  testa  using  con¬ 
vergence  criteria  aa  low  at  0.05X  ware  undertaken  for 
the  normalised  aum  of  residual  naas  source  magnitudes, 
fnr  example. 

Streamline  Plots 

Figures  3(a)  and  (b)  show  predicted  dimensionless 
streamline  patterns  calculated  and  plotted  by  coaputer . 
Tha  recirculation  call  given  la  Fig.  9(a)  la  tha  ex¬ 
pected  result  of  a  free  shear  layer  flowing  over  a 
cavity.  Tha  pradintiona  give  no  indication  that  an 
additional  tall  recirculation  a all  exists  in  tha 
corner  at  tha  basa  of  aaoh  tooth.  For  tha  case  of 
negligible  inlet  (l.a.  leakage)  flow  rata,  a  pair  of 
countar-rotating  Taylor-Ilka  eella  ara  axpaotad  within 
tha  cavity.  A  predicted  velocity  vector  plot  of  this 
flowfleld  is  presented  by  Stoff.  Ha  also  shows  a 
similar  plot  for  tha  seas  of  mall  inlet  flow  rata 
U/V  “  0.06.  That  case  exhibits  both  rselreulation 
calls,  tha  upstream  ona  being  extended  downatrean  near 
tha  free  shear  layer,  and  tha  other  occupying  only  tha 
lower  right  .hand  corner  of  tha  rectangular  cavity. 

Also  observe  in  Fig.  3(a)  that  the  dividing 
streamline  exhibits  a  r settlement  stagnation  point 
slightly  below  the  peripheral  corner  of  the  downstream 
tooth.  This  la  also  trua  of  tha  typical  cavity  design 
aa  seen  la  Fig.  3(b).  Tha  predicted  rselreulation 
strength  for  tha  typical  design  is  approxlnatsly  292 
of  tha  leakage  flow  rata,  wharaaa  it  la  39X  for  tha 
generic  case. 

Profile  Plots 

Many  quantitative  details  for  tha  flow  field  in  a 
labyrinth  cavity  ara  not  found  elsewhere,  particularly 
concerning  tha  inlet  boundary  values.  Therefore 
numerous  radial  profiles  of  aaoh  flow  variable  ara 
included  hare.  For  velocity  predictions,  three 
■apavate  solutions  ara  ahown,  employing i  (a)  tha 
QUICK  formulation  with  a  33  x  31  grid,  (b)  tha  Hybrid 
formulation  with  a  53  x  33  grid,  and  (c)  tha  Hybrid 
formulation  with  tha  33  x  31  grid.  Both  gride  are 
uniform  for  comparison  purposes. 

All  tha  data  presented  hare  was  calculated  and 
plotted  by  computer.  Figures  6,  8,  and  9  exhibit  tha 
distribution  of  predicted  axial,  radial,  and  awirl 
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(a)  Generic  cavity 
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(4)  Generic  cavity 
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(b)  Typical  cavity 

Fig.  3  Fredloted  dlaeasionlass  strasnlln*  pattern 
fori  (a)  a  lanarlo  and  (b)  a  typical 
labyrinth  aaal  cavity 


velocity  valuaa.  wharaaa  Fig.  7  provldaa  a  Magnified 
view  of  a  portion  of  Fig.  6  for  detailed  iuapeotiou. 
Figure  6(a)  above  the  dlaanalouleaa  Mean  axial 
velocity  distribution  for  the  generic  cavity.  The 
axial  velocity  value*  of  the  three  eolutlona  given  ere 
in  agreement.  Thee*  eolutlona  indicate  that  a  vary 
large  value  of  Su/tr  occurs  near  the  separation  cor¬ 
ner  j  this  la  Indicative  of  high  turbulenca  energy 
generation  In  the  free  shear  layer. 

The  magnified  view  (with  expanded  velocity  scale) 
of  Fig.  7(a)  reveals  thee,  Indeed,  the  33  x  31  QUICK 
and  33  x  33  Hybrid  solutions  are  saeentlally  Identical. 
This  Indicates  thee  both  solutions  era  accurate  and 
that  the  QUICK  fonulatlon  has  bean  correctly  inpls- 
nanted.  The  33  x  31  Hybrid  and  33  x  S3  Hybrid  aoiu- 
:,lone  differ  very  allghely. 

The  axial  velocity  distribution  for  the  typical 
cavity  In  Fig.  6(b'  la  very  similar  to  that  in  Fig. 
6(a).  The  neat  sigalf leant  differences  appear  at 
x/L  «  0.33,  both  In  the  free  shear  layer  end  in  the 
lower  axtranlty  of  the  cavity,  la  with  the  generic 
cavity,  the  Magnified  view  in  Fig.  7(b)  clsarJy  above 


(b)  Typical  cavity 


Fig.  6  Predicted  dlaeneioalese  naan  axial  velocity 
for  (a)  e  generic  and  (b)  a  typical  labyrinth 
seal  cavity  showing  33  x  31  QUICK  (———), 

33  x  33  Hybrid  < - ),  and  33  x  31 

Hybrid  (—  1  "  —  ■■■-»-■)  solutions 

that  tha  33  x  31  QUICK  and  the  33  x  S3  Hybrid  solu¬ 
tions  of  anlal  velocity  are  nearly  Identical.  Hota 
that  la  this  case  this  33  x  31  Hybtld  solution  is  not 
grid  independent  in  tha  free  shear  layer.  Tha  axial 
velocity  values  of  this  solution  devlats  froa  those  of 
tha  other  two  solutions  by  ss  Much  as  202  of  U. 

Observe  st  this  point  that  tha  QUICK  approach  - 
aa  inplanented  hart  -  sxhlblts  a  Major  acono«iaal 
advantage  over  tha  Hybrid  approach.  In  obtaining  a 
grid  independent  solution,  utilisation  of  the  Hybrid 
nethod  requires  S3  x  33  grid  points  (contusing  1.79  x 
10*  CPU  seconds  on  s  Print  730  coeputar),  wharaaa 
utilisation  of  QUICK  requires  only  33  x  31  points 
(using  only  7.79  x  103  CFU  seconds),  This  execution 
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Up.  *  Predicted  dlaanalonlaaa  maan  radial  velocity 
fort  (a)  a  generic  and  (b)  a  typical 
labyrinth  aaal  cavity  ihovini  33  x  31  QUICK 

( - ),  33  X  33  Hybrid  ( - - )  and 

33  x  31  hybrid  ( - —  - . — -)  solution* 


dacreaae  at  aanatrictlona  In  flow  area.  The  variation 
of  flow  area  la  Indicated  by  the  dividing  atraaallna 
In  Figure  3(a).  A  eharp  poaltive  preaaure  peak,  which 
occur*  In  an  annular  fashion  near  x/L  •  0.U3,  result* 
frow  the  flow  atagnatlon  on  the  downstream  tooth.  Tha 
resulting  overall  praiaure  drop  free  inlet  to  outlet 
of  thla  cavity  at  tha  stated  conditions  la  predicted 
ae  approximately  (.Ox  10*  kJa.  Only  the  33  x  31 
qUICK  solution  1*  presented  hare,  ae  value*  true  tha 
other  two  eolations  show  insignificant  deviations. 

The  preaeur*  distribution  within  the  typical 
cavity  1*  seen  in  Pig.  10(b)  tc  exhlbl  similar 
trend*.  However,  leakage  flow  pveature  varies  in  tha 
atveMMis*  direction  such  acre  gradually  for  thla 
gaoaetry  until  about  x/L  ■  0.(3.  Prooaading  down- 
■ersaa,  praaaura  dacraaaaj  eora  rapidly  than  for  tha 
generic  case  ae  tha  flow  arts  constrict*  la  passing 
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Pig.  9  Predicted  d learn lonlas*  eean  swirl  velocity 
fori  (a)  a  generic  end  (b)  a  typical 
labyrinth  aaal  cavity  showing  33  x  31  QUICK 

(  .  ) ,  33  x  33  Hybrid  (■-"  —  )  and 

33  x  31  Hybrid  ( - - )  solutions 


ever  tha  dovnitrsaa  tooth,  Tha  predicted  bulk 
praaaura  stop  for  thla  cavity  la  approximately  s.o  x 
10*  kPa,  an  incresaa  of  331  from  tha  i,enerlc  geonatry, 
(tease,  epf.tt  free  rotordynamic  stablllhj  coueidera- 
tiona,  the  typical  design  le  therefore  preferred  over 
its  generic  counterpart. 

As  nantionad  earlier,  tha  larg*  value  of  »u/#r  in 
th*  fro*  shear  layer  yields  a  high  gaaaratlon  rata  of 
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n«.  10  Predicted  dinenaicmleaa  Man  ralativa 
praaaura  ton  (a)  a  generic  and  (b)  a 
typical  labyrinth  anal  cavity  showing  tha 
33  x  31  QUICK  aolutlon 


turbulence  kinetic  energy.  this  of  couraa  promotes 
tha  large  ovatall  praaaura  drop.  Ilauoa,  tha  station 
naalaua  ef  turbulaoca  anargy  occurs  Mar  tha  traa 
abaav  ration  aa  saan  In  rigs.  11(a)  and  (b).  Hots 
that  tha  typical  cavity  gaoMtry  exhibits  a  signifi¬ 
cantly  higher  laval  of  turbulaoca  energy  than  for  tha 
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rig.  11  Predicted  uleaarlotUeaa  turbulence  kinetic 
energy  fori  (c)  a  generic  end  (b)  a 
typical  labyrinth  seal  cavity  shoving  tha 
33  x  31  QUICK  aolutlon 


generic  cast.  The  greatest  value  occurs  Is  both  cases 
neer  the  rssttschneat  stagnation  point.  1  significant 
portion  of  this  turbulence  energy  Is  convactad  over  the 
downatroa*  tootu  sa  indicated  by  tha  distrlbutloc  at 
x/L  *1.0.  la  expected  for  both  cavities  considered 
liars,  the  recirculation  region  la  assn  to  be  an 
sff active  turbulence  energy  sink. 

concurs  DIO  RUUJJT4 


Tha  coaplste  distrlbutloc  of  eauh  quantity 
revealing  tha  flov  character  vlthln  e  realistic  leby- 
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ri|.  7  Magnified  view  of  predicted  dimensionleat  Mas  axial  velocity  for 

(a)  a  generic  and  (b)  a  typical  labyrinth  aaal  cavity  a bovine  33  x  31 

QUICK  ( - ) ,  S3  x  S3  Hybrid  ( - )  and  33  x  31  Hybrid 

— -  solutions 


time  for  the  praaant  QUICK  varaion  of  tha  coda  la  only 
AA*  of  that  required  by  tha  Hybrid  varaion.  TUia  la 
•oatly  attributed  to  the  vail  known  apparent  sluggish 
raiponaa  to  (rid  refinement  of  tha  upwind  portion  of 
the  Hybrid  achaaa.  Aa  daw  narrated  by  Han  at  al.  [22], 
thla  lack  of  raaponaa  la  actually  a  talaa  diffualon 
affect. 

Aa  axpaotad.  there  la  nenllalbla  naan  radial 
velocity  In  the  leakage  flow  region  of  tha  ganarlo 
cavity  aa  aaan  in  Piga,  8(a)  and  (b)  axoapt  naar 
tha  atagnation  point  on  tha  downstream  tooth,  leal 
daalgnara  can  perhaps  reduce  tibia  leakag*  contribution 
by  altering  tha  tooth  ahapa  at  thla  location,  liaan- 
t tally  tha  asms  valuaa  raault  fro*  tha  three  aolutiona. 
Tha  avirl  velocity  in  rigs.  9(a)  and  (b)  la  nerval laid 
by  the  abaft  clraumfarentlal  velocity.  Tha  valuaa 


for  both  aaal  designs  are  almoat  unlforn  throughout 
tha  antira  doaaln.  They  era  0.65  and  0.69  for  Ugn. 
9(a)  and  (b)  respectively.  Noe  ahowc  ia  a  vary 
•harp  gradient  naar  each  surface.  All  three  aolutiona 
give  aaaantlally  the  ea*a  avirl  valoeley  valuaa  for 
tha  ganarlo  daaign  vhlla  the  33  x  31  Hybrid  solution 
la  not  quite  grid  independent  for  tho  typical  oaaa. 

Figures  10(a)  and  (b)  show  the  distribution  of 
praaaura  relative  to  tbs  cavity  Inlet  stator  vail 
pressure  ?ow>  In  tha  laakaga  flow  region  of  tha 
ganarlo  cavity,  praaaura  dacraaaaa  la  tha  itraananlaa 
direction  until  about  x/L  ■  0.3.  Than  it  lneraasea 
slightly  until  approximately  x/L  ■  0.78,  whereupon 
it  suddenly  dacraaaaa  aa  tha  flew  accelerate*  over  tha 
downstream  tooth.  Thla  pressure  variation  la  a 
reflection  of  tha  gemoulli  efface,  of  praaaura 
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rinth  seal  cevity  la  documented  hara  for  tha  firet 
time.  Thia  aarvaa  to  provide  aaal  designers  with  new 
inflight  regarding  important  flowfield  datalla.  Thaaa 
raaulta  further  allow  ona  to  advene*  tha  development 
of  simple  aapiclcally-baaad  algebraic  design  aodala  In 
currant  use.  Ona  raault  of  particular  Intaraat  In 
thla  ragard  la  that,  for  similar  caTltlaa ,  tha  diman- 
alonlaaa  awirl  velocity  may  ba  tccurataly  approximated 
by  a  unlfomly-dlatrlbutad  valua  of  0.69.  Anothar 
Important  finding  la  tha  shapa  of  tha  dividing  etre.ea- 
llna  and  how  It  Inf luancaa  laakaga  flow  praaaura 
variation. 

Tha  maaarieal  aapact  of  thla  Isvaatlgatlon  con'* 
atltutaa  a  tlaaly  and  atriagant  taat  of  tha  lnharant 
atabllity  of  QUICK  for  a  raaliatle  application.  Tha 
flowa  oonaldarad  hara  ara  highly  nonllnaar.  aa  the 
praaanca  of  awirl  momentum  la  known  [26]  to  laeraaoa 
tha  praaauro-valoclty  coupling  conaidarably.  Aa 
Inplanantad  hara,  QUICK  contributed  no  notleaabla 
atabllity  difficulty  for  aaeh  of  aavaral  aueh  highly 
nonllnaar  turbulent  flow  problana.  Moreover,  for  a 
given  laval  of  aesuraoy,  thla  varalon  of  tha  coda 
consumed  only  44X  of  tha  execution  tlna  required  by 
tha  Hybrid  varalon. 
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